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ABSTRACT 

Changes  in  Measured  Lightning  Return  Stroke  Peak  Current 
After  the  1994  National  Lightning  Detection  Network  Upgrade. 

(May  1997) 

Robert  Scott  Wacker,  B.S.,  United  States  Air  Force  Academy 
Chair  of  Advisory  Committee:  Dr.  Richard  E.  Orville 

Since  a  comprehensive  upgrade  of  the  US  National 
Lightning  Detection  Network  (NLDN)  in  1994,  the  mean  peak 
current  of  detected  cloud-to-ground  (CG)  lightning  flashes 
has  decreased,  the  number  of  detected  flashes  has  increased, 
and  the  percentage  of  flashes  with  positive  polarity  has 
increased.  NLDN  archived  data  for  the  years  1989-95  was 
analyzed  to  characterize  this  change.  Models  of  lightning 
return  stroke  radiation  and  its  detection  by  NLDN  sensors 
were  developed  to  simulate  the  effects  of  the  changes  made 
to  the  sensors  and  explain  the  cause  of  the  post-upgrade 
changes  in  NLDN  data. 

The  US  negative  mean  peak  current  decreased  from  a  pre¬ 
upgrade  (1989-93)  mean  of  37.5  kA  to  a  1995  value  of  30.2 
kA,  a  decrease  of  3.39  standard  deviations.  The  positive 
mean  peak  current  decreased  from  54.4  kA  to  31.6  kA,  a  4.97 
standard  deviation  decrease.  The  US  negative  flash  count 
increased  1.16  standard  deviations,  from  a  pre-upgrade  mean 


IV 


of  16.7  million  flashes  yr"^  to  20.6  million  flashes  in  1995. 
The  positive  flash  count  increased  6.22  standard  deviations, 
from  an  average  of  696  000  flashes  yr”^  before  the  upgrade  to 
2.1  million  flashes  in  1995.  Both  the  negative  and  positive 
flash  count  increases  were  predominant  at  low  peak  currents. 

A  model  of  return  stroke  detection  by  NLDN  sensors  was 
used  to  simulate  one  of  the  adjustments  made  to  NLDN  sensors 
during  the  upgrade.  Decreasing  the  pulsewidth  detection 
criterion  used  by  the  sensors  increases  their  effective 
detection  range,  which  increases  their  sensitivity  to  weak 
flashes  (due  to  NLDN  sensor  geometry,  increasing  sensitivity 
has  little  effect  on  detection  of  strong  flashes) .  The 
increased  detection  of  weak  flashes  accounts  for  the 
decrease  in  mean  peak  currents  and  the  increase  in  flash 


counts . 


CHANGES  IN  MEASURED  LIGHTNING  RETURN  STROKE  PEAK  CURRENT 
AFTER  THE  1994  NATIONAL  LIGHTNING  DETECTION  NETWORK  UPGRADE 

A  Thesis 
by 

ROBERT  SCOTT  WACKER 

Submitted  to  the  Office  of  Graduate  Studies  of 
Texas  A&M  University 

in  partial  fulfillment  of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 


May  1997 


Major  Subject:  Meteorology 


CHANGES  IN  MEASURED  LIGHTNING  RETURN  STROKE  PEAK  CURRENT 
AFTER  THE  1994  NATIONAL  LIGHTNING  DETECTION  NETWORK  UPGRADE 

A  Thesis 
by 

ROBERT  SCOTT  WACKER 

Submitted  to  Texas  A&M  University 
in  partial  fulfillment  of  the  requirements 
for  the  degree  of 

MASTER  OF  SCIENCE 


Approved  as  to  style  and  content  by: 


(Member)  (Head  of  Department) 


May  1997 

Major  Subject:  Meteorology 


Ill 


ABSTRACT 

Changes  in  Measured  Lightning  Return  Stroke  Peak  Current 
After  the  1994  National  Lightning  Detection  Network  Upgrade. 

(May  1997) 

Robert  Scott  Wacker,  B.S.,  United  States  Air  Force  Academy 
Chair  of  Advisory  Committee:  Dr.  Richard  E.  Orville 

Since  a  comprehensive  upgrade  of  the  US  National 
Lightning  Detection  Network  (NLDN)  in  1994,  the  mean  peak 
current  of  detected  cloud-to-ground  (CG)  lightning  flashes 
has  decreased,  the  number  of  detected  flashes  has  increased, 
and  the  percentage  of  flashes  with  positive  polarity  has 
increased.  NLDN  archived  data  for  the  years  1989-95  was 
analyzed  to  characterize  this  change.  Models  of  lightning 
return  stroke  radiation  and  its  detection  by  NLDN  sensors 
were  developed  to  simulate  the  effects  of  the  changes  made 
to  the  sensors  and  explain  the  cause  of  the  post-upgrade 
changes  in  NLDN  data. 

The  US  negative  mean  peak  current  decreased  from  a  pre¬ 
upgrade  (1989-93)  mean  of  37.5  kA  to  a  1995  value  of  30.2 
kA,  a  decrease  of  3.39  standard  deviations.  The  positive 
mean  peak  current  decreased  from  54.4  kA  to  31.6  kA,  a  4.97 
standard  deviation  decrease.  The  US  negative  flash  count 
increased  1.16  standard  deviations,  from  a  pre-upgrade  mean 


IV 


of  16.7  million  flashes  yr“^  to  20.6  million  flashes  in  1995. 
The  positive  flash  count  increased  6.22  standard  deviations, 
from  an  average  of  696  000  flashes  yr”^  before  the  upgrade  to 
2.1  million  flashes  in  1995.  Both  the  negative  and  positive 
flash  count  increases  were  predominant  at  low  peak  currents. 

A  model  of  return  stroke  detection  by  NLDN  sensors  was 
used  to  simulate  one  of  the  adjustments  made  to  NLDN  sensors 
during  the  upgrade.  Decreasing  the  pulsewidth  detection 
criterion  used  by  the  sensors  increases  their  effective 
detection  range,  which  increases  their  sensitivity  to  weak 
flashes  (due  to  NLDN  sensor  geometry,  increasing  sensitivity 
has  little  effect  on  detection  of  strong  flashes) .  The 
increased  detection  of  weak  flashes  accounts  for  the 
decrease  in  mean  peak  currents  and  the  increase  in  flash 


counts . 


V 


DEDICATION 

This  thesis  is  dedicated  to  my  parents.  Herb  and  Judy 
Wacker,  for  the  lifetime  of  love,  encouragement,  and  support 
they  have  given  me. 


VI 


ACKNOWLEDGMENTS 

The  author  wishes  to  acknowledge  the  Fannie  and  John 
Hertz  Foundation  for  its  financial  support  of  his  graduate 
studies;  the  United  States  Air  Force  for  providing  the 
opportunity  to  pursue  his  master's  degree;  and  Dr.  Richard 
E.  Orville  and  the  Cooperative  Institute  for  Applied 
Meteorological  Studies  at  Texas  A&M  University  for  the 
guidance  and  resources  necessary  to  complete  it. 


vii 


TABLE  OF  CONTENTS 

Page 

ABSTRACT .  iii 

DEDICATION .  v 

ACKNOWLEDGMENTS .  vi 

TABLE  OF  CONTENTS .  vii 

LIST  OF  FIGURES . . ix 

LIST  OF  TABLES .  xii 

CHAPTER 

I  INTRODUCTION .  1 

II  BACKGROUND .  3 

1.  The  lightning  return  stroke .  4 

2.  Wideband  magnetic  direction  finder 

operation .  12 

3.  The  National  Lightning  Detection 

Network .  16 

III  MEASURING  PEAK  CURRENT .  22 

1.  Electromagnetic  theory .  22 

2.  Peak  current  calibration .  31 

IV  NLDN  DATA  ANALYSIS .  35 

1.  Previous  work .  35 

2.  Data  set  and  analysis  methods .  35 

3.  Results .  38 

V  RETURN  STROKE  DETECTION  MODELING .  71 

1.  Previous  work .  71 

2.  Finite  difference  return  stroke 

field  model . .  7  6 

3.  Stochastic  return  stroke  detection 

model .  88 

VI  DISCUSSION  AND  CONCLUSIONS .  96 


viii 


REFERENCES . 
APPENDIX  A: 

APPENDIX  B: 
APPENDIX  C: 
APPENDIX  D: 


1.  Statistical  significance  of  data...  96 

2.  Detection  modeling  results. .  97 

3.  Possible  explanation  of  peak 

current  decrease  and  flash  count 
increase . 98 

4.  Recommendations .  101 


.  103 

ORVILLE  (1991)  AND  IDONE  ET  AL.  (1993) 


PEAK  CURRENT  CALIBRATIONS .  106 

NLDN  DATA  ANALYSIS  SAMPLE . 119 

FINITE  DIFFERENCE  FIELD  MODEL .  127 

STOCHASTIC  DETECTION  MODEL .  132 


VITA 


139 


IX 


LIST  OF  FIGURES 

FIGURE  Page 

2.1  Chronology  of  a  CG  lightning  flash .  7 

2.2  Magnetic  field  loop  antenna  and  electric 

field  flat  plat  antenna  configuration  in 
wideband  magnetic  direction  finders .  13 

2.3  The  National  Lightning  Detection  Network 

after  the  1994  upgrade... .  19 

3.1  Typical  first  and  subsequent  return  stroke 

waveforms  with  dart-stepped  leaders  and  dart 
leaders .  29 

4.1  Kansas  annual  and  monthly  negative  flash 

count,  1989-95 .  41 

4.2  Kansas  annual  and  monthly  negative  mean  peak 

current,  1989-95 .  42 

4.3  Kansas  annual  and  monthly  positive  flash 

count,  1989-95 .  43 

4.4  Kansas  annual  and  monthly  positive  mean  peak 

current,  1989-95 .  44 

4.5  Ohio  Valley  annual  and  monthly  negative 

flash  count,  1989-95 .  47 

4.6  Ohio  Valley  annual  and  monthly  negative  mean 

peak  current,  1989-95 .  48 

4.7  Ohio  Valley  annual  and  monthly  positive 

flash  count,  1989-95 .  49 

4.8  Ohio  Valley  annual  and  monthly  positive  mean 

peak  current,  1989-95 .  50 

4.9  Florida  annual  and  monthly  negative  flash 

count,  1989-95 .  53 

4.10  Florida  annual  and  monthly  negative  mean 

peak  current,  1989-95 .  54 

4.11  Florida  annual  and  monthly  positive  flash 

count,  1989-95 .  55 


X 


FIGURE  Page 

4.12  Florida  annual  and  monthly  positive  mean 

peak  current/  1989-95 .  56 

4.13  US  annual  and  monthly  negative  flash  count, 

1989-95 .  59 

4.14  US  annual  and  monthly  negative  mean  peak 

current,  1989-95 .  60 

4.15  US  annual  and  monthly  positive  flash  count, 

1989-95 .  61 

4.16  US  annual  and  monthly  positive  mean  peak 

current,  1989-95 .  62 

4.17  US  negative  peak  current  distribution, 

1989-95 .  65 

4.18  US  negative  pre-  to  post-upgrade  flash  count 

ratio .  66 

4.19  US  positive  peak  current  distribution, 

1989-95 .  67 

4.20  US  positive  pre-  to  post-upgrade  flash  count 

ratio .  68 

5.1  Finite  difference  field  model  channel 

geometry .  78 

5.2  Model  channel  current  vs.  height  15  p.s  after 

return  stroke  initiation .  80 

5.3  Electric  field  strength  and  magnetic  flux 

density  measured  100  km  from  the  30  kA  peak 
current  model  return  stroke  in  Fig.  5.2 .  82 

5.4  Raw  and  normalized  signal  strength  vs. 

range .  84 

5.5  Signal  pulsewidth  vs.  range  for  the  30  kA 

return  stroke  modeled  in  Fig.  5.2 .  85 

5.6  Detection  range  vs.  peak  current  for  modeled 
return  strokes  having  waveforms  identical  to 
Fig.  5.2  using  detection  pulsewidth  criteria 
of  10  |j,s  and  7^is 


87 


xi 


FIGURE  Page 

5.7  Detection  model  randomly  generated  return 

stroke  locations .  89 

5.8  Modeled  return  stroke  detected  flash  count 
vs .  peak  current  and  detected  flash  counts 

at  10  ps  and  7  ^s  pulsewidth  criteria .  92 


5.9  Model  detection  efficiency  vs.  peak  current 

at  10  {IS  and  7  jis  pulsewidth  criteria  and 
model  post-upgrade  flash  count  vs.  peak 
current  normalized  to  pre-upgrade  flash 
count . 


93 


Xll 


LIST  OF  TABLES 

TABLE  Page 

4.1  Data  analysis  geographical  subsets .  36 

4.2  Post-  vs.  pre-upgrade  negative  flash  count 

and  mean  peak  current  comparison .  38 

4.3  Post-  vs.  pre-upgrade  positive  flash  count 

and  mean  peak  current  comparison .  39 


1 


CHAPTER  I 
INTRODUCTION 

Since  1989  the  National  Lightning  Detection  Network 
(NLDN)  has  detected  cloud-to-ground  (CG)  lightning  flashes 
within  the  continental  United  States.  The  network  is 
administered  by  Global  Atmospherics,  Inc.,  of  Tucson, 
Arizona,  and  consists  of  approximately  one  hundred  lightning 
sensors  located  throughout  the  country.  These  sensors  are 
linked  via  satellite  to  a  centralized  network  control  center 
(NCC)  in  Tucson,  which  disseminates  detected  CG  flash  data 
to  NLDN  subscribers. 

NLDN  data  has  a  number  of  applications:  protection  of 
electrical  utility  assets;  forest  fire  prevention; 
protection  of  aircraft  and  spacecraft;  public  safety 
(according  to  Orville  1993,  lightning  is  the  largest  weather 
related  killer  in  the  United  States,  killing  100-150  people 
annually) ;  and  meteorological  research. 

During  the  latter  half  of  1994,  the  NLDN  underwent  a 
system-wide  upgrade  intended  to  improve  detection 
efficiency,  location  accuracy,  data  delivery,  and  network 
reliability.  This  upgrade  consisted  of  adjustments  to 
existing  sensor  hardware,  decommissioning  or  movement  of 
selected  sensors,  inclusion  of  a  new  type  of  sensor  absorbed 

This  thesis  follows  the  style  and  format  of  Monthly  Weather 
Review. 


from  another  network,  and  implementation  of  a  new  flash 
location  algorithm  (Cummins  et  al.  1995) . 
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Since  the  upgrade,  the  following  changes  have  been 
observed  in  NLDN  data:  1)  a  decrease  in  the  mean  peak 
current  of  detected  CG  flashes,  2)  an  increase  in  the  number 
of  detected  flashes,  and  3)  an  increase  in  the  percentage  of 
the  flashes  having  positive  polarity. 

The  purpose  of  this  thesis  is  to  describe 
quantitatively  the  post-upgrade  change  in  NLDN  peak  current 
and  flash  count  data  and  determine  its  probable  cause. 
Chapter  II  includes  background  on  CG  lightning  flashes, 
magnetic  direction  finder  (DF)  operation,  and  the  NLDN. 
Chapter  III  discusses  the  theoretical  basis  for  remote 
measurement  of  return  stroke  peak  current  and  the 
calibration  of  the  sensors  used  in  the  NLDN.  Chapter  IV 
contains  a  quantitative  analysis  of  the  changes  in  mean  peak 
currents  and  flash  counts  in  the  NLDN  archived  data. 

Chapter  V  presents  two  models  representing  the 
electromagnetic  radiation  emitted  by  a  CG  lightning  return 
stroke  and  its  detection  by  NLDN  sensors  before  and  after 
simulating  the  changes  made  to  them  during  the  NLDN  upgrade. 
Finally,  Chapter  VI  draws  conclusions  as  to  the  cause  of  the 
change  in  NLDN  data  based  on  the  similarities  between  the 
analysis  in  Chapter  IV  and  the  model  results  in  Chapter  V. 
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CHAPTER  II 
BACKGROUND 

The  purpose  of  this  chapter  is  to  provide  the  necessary 
background  knowledge  of  1)  the  physical  processes  involved 
in  lightning  flashes,  2)  the  operation  of  magnetic  direction 
finders,  and  3)  the  operation  of  the  NLDN  required  for 
understanding  of  the  later  chapters  discussing  radiation 
theory,  data  analysis,  and  modeling. 

1.  Terminology  and  coordinates 

Several  conventions  are  used  in  the  lightning 
literature  and  this  thesis.  The  first  distinguishes 
lightning  strokes  from  lightning  flashes.  A  stroke  is  a 
single  lightning  discharge.  A  lightning  flash  consists  of 
multiple  lightning  stokes  occurring  in  rapid  succession 
within  the  same  channel  from  cloud  to  ground.  Lightning 
flashes  appear  to  flicker  because  the  unaided  eye  is  just 
able  to  resolve  the  individual  luminosity  pulses  of  each 
stroke.  This  is  referred  to  as  the  flash's  multiplicity. 

Polarity  refers  to  the  sign  of  charge  a  flash  lowers  to 
ground.  Since  positive  and  negative  flashes  are  different 
in  several  respects,  they  are  considered  separately  in 
analyzing  lightning  data. 
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2 .  The  Lightning  Return  Stroke 

a.  Types  of  lightning  discharges 

Cloud-to-ground  (CG)  lightning  flashes  transfer  charge 
between  cloud  and  ground.  These  flashes  are  the  type  the 
NLDN  is  designed  to  detect  exclusively.  They  can  have  peak 
currents  of  tens  or  even  hundreds  of  kA. 

Cloud  flashes  occur  between  charge  regions  within  a 
cloud,  between  clouds,  or  between  a  cloud  and  the 
surrounding  air.  Cloud  flashes  comprise  approximately  80% 
of  all  lightning  discharges  and  have  peak  currents  of  only  a 
few  kA  (Orville  1993) .  The  electronics  of  the  sensors  used 
in  the  NLDN  are  specifically  designed  to  screen  out  cloud 
flashes . 

Triggered  flashes  occur  between  cloud  charge  regions 
and  man-made  objects  on  the  ground  or  in  the  air  which 
artificially  intensify  the  local  electric  field.  The  sensor 
calibration  discussed  in  Chapter  III  was  accomplished  using 
lightning  flashes  triggered  by  rockets  launched  into 
thunderstorms  at  the  Kennedy  Space  Center. 

Niomerous  other  names  exist  to  describe  lightning  or 
lightning-like  discharges;  heat  lightning,  bead  lightning, 
ribbon  lightning,  and  ball  lightning,  for  instance.  Some  of 
these  phenomena,  like  ribbon  lightning,  are  well  understood, 
while  others,  like  ball  lightning,  are  not.  Recently,  faint 


transient  optical  emissions  called  sprites  and  jets  have 
been  discovered  above  thunderstorms. 
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b.  CG  lightning  flash  processes 

A  CG  lightning  flash  occurs  at  the  end  of  a  complex 
chain  of  events  set  in  motion  by  the  nature  of  the  Earth- 
atmosphere  global  electrical  circuit. 

1)  The  Earth-atmosphere  circuit 

In  fair  weather,  the  ambient  electric  field  is  directed 
downward  with  a  strength  of  approximately  100  V  m~^  at  the 
earth's  surface,  decreasing  with  height.  The  total 
potential  difference  between  the  surface  and  50  km  is 
approximately  3  x  10®  V.  This  potential  difference  is  a 
consequence  of  the  net  negative  charge  of  approximately  4.5 
X  10®  C  residing  on  the  Earth's  surface.  The  Earth  and 
atmosphere  act  as  a  "leaky  spherical  capacitor",  with  a 
total  current  of  1000-1500  A  continuously  flowing  toward  the 
surface. 

Radioactive  decay,  cosmic  radiation,  and  short-wave 
ultraviolet  radiation  all  create  charge  in  the  atmosphere 
through  ionization  (Uman  1987) .  However,  we  would  expect 
the  downward  current  to  neutralize  the  surface  charge  and 
eliminate  the  fair  weather  electric  field  in  a  matter  of 
minutes.  Therefore,  there  must  be  some  mechanism  which 
continually  maintains  the  ambient  charge  distribution. 
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Thunderstorms  serve  this  role  by  carrying  negative  charge 
from  cloud  to  ground,  thus  acting  as  a  battery  in  the  Earth- 
atmosphere  circuit. 

2)  Charge  separation 

Fig.  2.1  illustrates  the  processes  occurring  before  and 
during  a  lightning  flash.  The  first  of  these  is  charge 
separation  within  a  cloud  (t=0  in  Fig.  2.1).  The  exact 
mechanism  of  charge  separation  is  not  well  understood,  but 
it  is  thought  to  be  a  consequence  of  collisions  between 
mixed-phase  precipitation  (i.e.  rain  and  graupel)  within  the 
updraft.  Due  to  the  requirement  for  both  frozen  and  liquid 
precipitation,  the  level  where  charge  separation  occurs  is 
dependent  upon  the  height  of  the  freezing  level. 

A  typical  cumulonimbus  cloud  contains  a  double-dipole 
charge  structure:  a  large  region  of  positive  charge, 
perhaps  40  C,  exists  near  the  -45°C  level;  another  large 

region  of  equal  negative  charge  exists  lower,  near  the  -15°C 
level;  and  a  small  region  of  positive  charge  typically 
exists  near  cloud  base  (Uman  1987) .  This  charge 
distribution  induces  a  positive  charge  region  at  the  surface 
below  the  cloud  and  results  in  a  potential  difference 
between  cloud  base  and  ground  on  the  order  of  10^  V.  Within 
this  large-scale  structure,  smaller  sub-regions  of  cloud 
charge  exist. 


Fig.  2.1.  Chronology  of  a  CG  lightning  flash  (adapted  from 
Uman  1987) .  Note  the  approximate  time  reference  below  each 
frame. 
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3)  Stepped  leader 

When  the  electric  field  between  cloud  charge  regions 
exceeds  the  ionization  potential  of  air  (approximately  500 
000  V  depending  on  pressure),  breakdown  occurs  (t=1.00 
ms  in  Fig.  2.1).  The  ionized  air  forms  a  conducting  channel 
into  which  cloud  charge  flows,  further  intensifying  the 
potential  gradient,  and  hence  the  electric  field.  This 
causes  another  breakdown  to  occur,  adding  another  ionized 
segment  to  the  conducting  channel.  This  faintly  luminous 
process,  the  stepped  leader,  repeats  itself  approximately 
every  50  |xs,  adding  segments  tens  of  m  long  to  the  channel 
and  branching  occasionally  as  it  moves  toward  ground  with  an 
average  velocity  of  approximately  2  x  10®  m  s”^  {t=1.00  ms  to 
t=19.00  ms  in  Fig.  2.1). 

4)  Return  stroke 

As  the  stepped  leader  nears  ground,  the  potential 
difference  between  the  tips  of  the  leader  branches  and 
ground  remains  constant,  but  as  the  distance  decreases  the 
potential  gradient  grows  very  large,  especially  in  the 
presence  of  tall,  pointed  objects  on  the  ground  where  charge 
can  accumulate.  Breakdown  occurs  at  such  an  object  when  the 
field  grows  large  enough,  inducing  an  upward  streamer  (Fig. 
2.1,  t=20.00  ms).  When  the  upward  streamer  meets  the  tip  of 
one  of  the  stepped  leader  branches  the  cloud  charge  is 
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short-circuited  to  ground,  initiating  the  return  stroke  (Fig 
2.1,  t=20.10  ms  to  20.20  ms).  A  wave  of  ground  potential 
travels  up  the  channel  toward  the  cloud  charge  region  at  a 
velocity  of  approximately  10®  m  s”^.  Once  this  potential 
wave  passes  a  point  on  the  channel,  the  charge  stored  there 
flows  to  ground.  Charge  stored  in  a  branch  flows  to  ground 
after  the  potential  wave  passes  the  point  where  it  connects 
to  the  channel.  The  resulting  large  current  flowing  in  the 
channel  ionizes  the  various  gaseous  constituents  of  the 
atmosphere,  producing  the  intense  luminosity  we  see.  It 
also  heats  the  channel  to  approximately  30  000  K,  causing  an 
overpressure  of  about  100  atm.  This  causes  the  outward- 
propagating  shockwave  heard  as  thunder. 

5)  Dart  leader 

The  first  return  stroke  discharges  one  region  of  charge 
within  the  cloud.  The  channel,  however,  remains  a  conductor 
until  the  ionized  air  composing  it  can  recombine.  While 
this  is  taking  place,  the  channel  may  connect  to  another 
cloud  charge  region  (Fig.  2.1,  t=40.00  ms).  A  new  leader, 
the  dart  leader,  moves  down  the  partially-recombined 
channel,  recharging  it.  Since  it  already  has  the  existing 
channel  to  follow,  the  dart  leader  about  ten  times  faster 
than  the  stepped  leader  (Fig.  2.1,  t=60.00  ms  to  62.00  ms). 
Sometimes  recombination  will  have  already  occurred  at  the 
bottom  of  the  channel  where  pressure  is  higher,  causing  the 
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dart  leader  to  finish  its  downward  movement  as  a  stepped 
leader. 

6)  Subsequent  return  stroke (s) 

Once  the  dart  leader  has  recharged  the  channel  and 
contacted  ground  another  return  stroke  occurs  (Fig.  2.1, 
62.05  ms) .  Subsequent  return  strokes  differ  from  first 
return  strokes  in  that  they  have  no  branches,  since  only  the 
existing  main  channel  is  re-ionized  by  the  dart  leader.  The 
dart  leader-subsequent  return  stroke  process  can  repeat 
itself  several  times  with  tens  of  ms  between  return  strokes 
(Uman  1987) .  In  a  significant  fraction  of  CG  return 
strokes,  a  continuing  current  of  approximately  100  A  may 
flow  for  as  long  as  tenths  of  a  second  (Uman  1987) . 

7)  Return  stroke  characteristics 

Typical  return  stroke  characteristics  are  summarized 

here : 

Time  required  for  the  return  stroke  to  propagate  up  the 
channel  can  range  from  10  ^is  (Beasley  et  al.  1993)  to  20-50 
}is  (Cummins  et  al.  1995)  to  100  |4,s  (Uman  1987).  Time 
between  subsequent  return  strokes  ranges  between  20  and  100 
ms  (Cummins  et  al.  1995). 

Return  stroke  velocity  averages  10®  m  s‘^  (Beasley  et 
al.  1995)  and  decreases  as  each  branch  is  passed  (Uman 


1987)  . 
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Cummins  et  al.  (1995)  report  a  peak  current  range  of  5 
to  300  kA.  Mean  peak  current  appears  to  decrease  with 
latitude  (Orville  1990) . 

Multiplicity  averages  between  2  and  3,  but  values  may 
range  between  1  and  20  (Cummins  et  al.  1995) . 

Orville  (1994)  reports  a  total  CG  flash  count  of 
approximately  10^  CG  flashes  per  year  in  the  United  States. 
About  4%  of  CG  flashes  are  of  positive  polarity,  with 
considerable  latitudinal  and  seasonal  variation, 
c.  Positive  flashes 

CG  flashes  may  lower  positive  or  negative  charge  to 
ground.  Positive  CG  flashes,  which  comprise  approximately 
5%  of  all  CG  flashes,  differ  significantly  from  negative 
flashes.  Typically  positive  flashes  have  higher  peak 
currents  (approximately  45  )cA  mean  peak  current  vs. 
approximately  30  kA  for  negative  flashes)  and  transfer  more 
charge.  Positive  leaders  typically  are  not  branched  and 
positive  flashes  usually  do  not  have  subsequent  return 
strokes . 

The  percentage  of  positive  flashes  is  higher  in  winter, 
at  high  latitudes,  at  high  elevations,  and  in  severe  storms 
(Uman  1987;  Orville  1994) .  This  leads  us  to  conclude  that 
positive  CG  flashes  occur  between  the  upper  (positive) 
charge  region  of  the  cloud  and  ground.  A  lower  freezing 
level  or  higher  terrain  decrease  the  distance  from  the 
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positive  charge  to  ground,  and  strong  shear  tilts  the  dipole 
structure,  providing  a  more  direct  path  from  the  positive 
charge  to  ground. 

3.  Wideband  magnetic  direction  finder  op)eration 

An  entire  CG  lightning  flashes  can  last  as  long  as  1  s, 
while  the  individual  processes  within  it  occur  on  a  sub¬ 
microsecond  timescale,  so  the  frequencies  of  radiation 
emitted  during  a  flash  can  range  from  below  1  Hz  to  well  in 
excess  of  10  MHz  (Uman  1987) .  Therefore,  sensors  designed 
to  detect  lightning  radiation  waveforms  must  operate  over  a 
wide  frequency  range . 

The  first  wideband  magnetic  direction  finders  were 
developed  in  the  late  1970s  for  use  by  the  US  Bureau  of  Land 
Management  (BLM)  for  forest  fire  protection  in  the  Western 
United  States  (Krider  et  al.  1980).  The  descendants  of 
those  original  DFs  are  in  service  today  in  the  NLDN. 

For  reasons  discussed  in  Chapter  III,  these  DFs  are 
designed  to  measure  the  vertical  component  of  the  local 
electric  field  using  horizontal  flat  plate  antennae  and  the 
horizontal  components  of  the  local  magnetic  flux  density 
using  orthogonal  loop  antennae  in  the  vertical  plane.  Fig. 
2.2  depicts  the  antenna  configuration. 


Fig.  2.2.  Magnetic  field  loop  antenna  and  electric  field 
flat  plat  antenna  configuration  in  wideband  magnetic 
direction  finders  (top,  adapted  from  Lightning  Location  and 
Protection,  Inc.,  1992).  An  IMPACT  DF  manufactured  by  LLP, 
Inc.,  used  in  the  TOGA  COARE  field  experiment  (bottom). 
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a.  Magnetic  field  loop  antennae 

From  the  integral  form  of  Faraday's  law,  e  =  -dO/dt,  we 
see  that  a  changing  magnetic  flux,  O,  through  a  current  loop 

produces  an  electromotive  force,  e,  in  the  loop  (Griffiths 
1989) .  Two  orthogonal  vertical  loops  can  then  measure  the 
component  of  the  horizontal  magnetic  flux  normal  to  each 
loop  by  integrating  in  time  the  voltage  in  the  loop  (Krider 
et  al.  1980)  .  The  vector  Siam  of  the  two  components  yields 
the  total  horizontal  magnetic  flux  density.  Additionally, 
the  ratio  of  the  two  components  is  the  tangent  of  the 
azimuth  to  the  radiation  source  (Uman  1987) . 

b.  Electric  field  flat  plate  antenna 

Since  reciprocal  angles  have  the  same  tangent,  the 
vertical  electric  field  must  be  measured  to  determine  the 
polarity  of  the  return  stroke  and  eliminate  the  ambiguity  in 
the  azimuth  determined  by  the  vertical  loop  antennae.  This 
is  accomplished  using  an  elevated  horizontal  flat  plate. 

From  the  integral  form  of  Gauss'  law,  Q  =  EnSoA,  the  charge, 

Q,  induced  on  the  plate  is  proportional  to  the  electric 
field  component  normal  to  it,  En  (Griffiths  1989) .  A  change 
in  normal  electric  field  strength,  dEn/dt,  will  result  in  a 
current,  dQ/dt,  flowing  to  or  from  the  plate,  which  is 
integrated  to  yield  the  total  field  (Uman  1987) . 
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c.  Detection  criteria 

The  electronics  of  a  DF  are  configured  to  distinguish 
the  unique  radiation  waveforms  of  CG  return  strokes  from  the 
radiation  of  cloud  flashes  and  background  noise.  The 
waveform  characteristics  used  to  distinguish  a  return  stroke 
signature  are:  1)  its  risetime  from  the  detection  threshold 
to  the  initial  peak  value,  2)  its  pulsewidth,  3)  the 
magnitude  of  any  subsidiary  peaks,  and  4)  the  magnitude  of 
the  field  overshoot  (Krider  et  al.  1980) .  Return  stroke 
radiation  waveforms  are  discussed  in  greater  detail  in 
Chapter  III. 

d.  LLP  units 

The  magnetic  DFs  used  in  the  NLDN  are  manufactured  by 
Lightning  Location  and  Protection,  Inc.  (LLP) .  Their  output 
for  both  magnetic  flux  density  and  electric  field  strength 
are  calibrated  in  LLP  units.  One  LLP  unit  corresponds  to  a 
magnetic  flux  density  of  1.5  x  10”^°  wb  m”^  and  an  electric 
field  strength  of  4.5  x  10'^  V  m"^  (Hiscox  and  Cummins  1993). 
All  signal  strength  data  used  by  the  NLDN  are  in  LLP  units 
and  in  Chapter  III  a  calibration  of  return  stroke  peak 
current  in  terms  of  range  normalized  signal  strength  (in  LLP 
units)  is  presented. 
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4 .  The  National  Lightning  Detection  Network 

a.  History 

The  National  Lightning  Detection  Network  has  its  roots 
in  three  smaller  networks  that  came  into  existence  during 
the  1970s  and  '80s.  The  first  was  the  BLM  network 
established  in  the  late  1970s  for  forest  fire  prevention  in 
the  Western  United  States  and  Alaska  and  developed  by  Krider 
et  al.  (1980) .  Another  network  was  operated  by  the  National 
Severe  Storms  Laboratory  (NSSL)  in  the  Midwestern  United 
States.  The  first  network  with  centralized  data  processing 
was  developed  at  the  State  University  of  New  York  at  Albany 
(SUNYA)  and  initially  consisted  of  ten  sensors  covering  an 
area  along  the  East  Coast  from  North  Carolina  to  Canada 
(Cummins  et  al.  1995;  Orville  et  al.  1983). 

In  1987  the  three  separate  networks  were  combined  into 
the  NLDN  and  by  1989  enough  sensors  had  been  added  to 
provide  complete  coverage  of  the  continental  United  States 
(Cummins  et  al.  1995) .  In  1991,  GeoMet  Data  Services  (GDS) 
was  established  to  administer  the  NLDN  as  a  commercial 
venture.  GDS,  LLP,  and  Atmospheric  Research  Systems  (ARS, 
another  sensor  manufacturer)  comprise  Global  Atmospherics, 
Inc.,  with  headquarters  in  Tucson,  Arizona. 

b.  Flash  location  methods 

Flashes  are  located  by  the  NLDN  using  three  different 


methods : 
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1)  Time-of-arrival  (TOA)  method 

The  TOA  method  is  used  by  LPATS  sensors  developed  by 
ARS  and  incorporated  into  the  NLDN  in  the  1994  upgrade. 

LPATS  sensors  only  record  the  arrival  time  of  the  peak 
signal  strength  of  the  detected  return  stroke  signal;  they 
have  no  direction-finding  capability.  The  difference  in 
TOAs  recorded  by  a  pair  of  sensors  is  used  to  compute  a 
hyperbola  along  which  the  flash  must  have  occurred.  The 
intersection  of  hyperbolae  computed  from  multiple  pairs  of 
sensors  defines  the  flash  location.  At  least  four  sensors 
are  required  to  obtain  an  unambiguous  flash  location 
(Orville  1993) . 

2)  Triangulation  method 

When  azimuth  data  from  at  least  two  DFs  is  available, 
flash  locations  can  be  determined  by  triangulation.  This 
gives  good  results  as  long  as  the  flash  is  not  near  the 
baseline  between  two  DFs  whose  azimuths  are  being  used 
(Orville  1993) .  In  this  case,  azimuth  error  results  in  very 
large  position  error. 

3)  IMPACT  (Improved  Performance  from  Combined 
Technology)  method 

To  accommodate  the  addition  of  TOA  sensors  to  the 
existing  network  of  magnetic  DFs,  Global  Atmospherics 
developed  a  flash  location  algorithm  using  arbitrary 
combinations  of  TOA  and  azimuth  data.  Such  a  combination 
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defeats  the  shortcomings  of  each  method  used  individually 
(Cummins  et  al.  1995;  Orville  1993) . 

When  multiple  sensors  detect  a  flash  the  redundant  data 
is  used  to  optimize  the  flash  location  estimation.  As  a 
minimum,  the  IMPACT  method  requires  one  more  data  point  than 
the  number  of  parameters  desired.  For  instance,  a  two- 
dimensional  flash  location  and  time  requires  four  data 
points  (i.e.  two  azimuths  and  two  TOAs,  one  azimuth  and 
three  TOAs,  etc.)  (Cummins  et  al.  1995). 
c.  Network  configuration 

Since  the  1994  upgrade  the  network  consists  of  45 
existing  magnetic  DFs  updated  to  the  IMPACT  standard,  which 
provide  azimuth  and  TOA  based  on  the  Global  Positioning 
System  (GPS)  time  reference,  and  60  LPATS  sensors  acquired 
from  the  ARS  network  and  added  to  the  network,  which  provide 
GPS-based  time  of  arrival  only.  Fig.  2.3  depicts  the 
locations  of  both  types  of  sensors. 

The  sensors  send  their  data  via  satellite  link  to  the 
Network  Control  Center  (NCC)  in  Tucson.  At  the  NCC,  data 
from  all  sensors  in  the  network  is  combined  to  determine 
flash  time,  location,  peak  current,  polarity,  and 
multiplicity.  The  processed  flash  data  is  then  disseminated 
to  NLDN  subscribers  via  satellite.  The  entire  process  takes 
30-40  s  (Cummins  et  al .  1995). 


Fig.  2.3.  The  National  Lightning  Detection  Network  after 
the  1994  upgrade.  Red  triangles  represent  IMPACT  DFs  and 
blue  circles  represent  LPATS  time-of-arrival  (TOA)  sensors 
(adapted  from  Cummins  et  al.  1995) . 
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d.  The  1994  NLDN  upgrade 

Since  Global  Atmospherics  is  a  commercial  venture, 
changes  and  upgrades  to  improve  network  performance  have 
been  ongoing  since  the  NLDN's  inception.  For  example,  in 
late  1992  four  IMPACT  DFs  were  installed  in  Tennessee  and 
during  March-April  1993  the  original  DFs  from  the  BLM 
network  were  replaced  by  IMPACT  sensors  while  others  were 
added  in  South  Carolina  (Cummins  1995) . 

The  1994  upgrade,  however,  was  far  greater  in  scope. 
According  to  a  Global  Atmospherics  technical  note  describing 
the  upgrade,  its  primary  objective  was  to  improve  location 
accuracy  through  implementation  of  the  IMPACT  method.  It 
also  sought  to  improve  the  data  processing  infrastructure  to 
deliver  stroke  and  flash  data  in  real-time,  improve 
detection  efficiency  of  weak  flashes  (as  low  as  5  kA)  ,  and 
improve  the  long-term  reliability  of  the  NLDN  hardware 
(Cummins  et  al.  1995).  It  took  place  during  July-December 
1994,  with  the  majority  of  the  work  being  accomplished 
during  July  and  August  (Cummins  1996) . 

In  keeping  with  the  goal  of  increasing  weak  flash 
detection  efficiency,  several  hardware  changes  were  made  to 
NLDN  DFs:  1)  sensor  gain  was  increased  50%,  improving 
nominal  detection  range  from  400  km  to  600  km  (Cummins 
1996);  2)  the  pulsewidth  criterion  for  recognizing  return 
strokes  was  reduced  from  10  ps  to  7  jis  (Lucas  and  Orville 
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1996;  Cummins  1996) ;  and  3)  the  subsidiary  peak  tolerance 
was  increased  from  .85  of  the  initial  signal  strength  to  1.2 
(Cummins  1996) ,  allowing  detection  of  distant  weak  flashes, 
whose  ground  wave  may  be  attenuated  to  a  strength  comparable 
to  the  reflected  sky  wave.  The  effects  of  these  detection 
criteria  changes  are  discussed  further  in  Chapter  V. 
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CHAPTER  III 

MEASURING  PEAK  CURRENT 

The  purpose  of  this  chapter  is  to  discuss  the  physical 
and  mathematical  basis  for  measuring  peak  current  with 
remote . lightning  sensors.  It  presents  the  mathematical 
relation  between  return  stroke  peak  current  and  peak 
measured  field  strength  and  the  assumptions  used  to  derive 
it.  Finally/  it  presents  results  of  previous  sensor  peak 
current  calibrations/  corrects  for  an  erroneous  data  point 
in  one  calibration/  and  investigates  the  effects  of  using  a 
more  physically  reasonable  calibration. 

1 .  Electromagnetic  theory 
a.  Maxwell's  equations 

Electromagnetic  radiation  is  governed  by  Maxwell's 
equations  (Equations  3.1a  through  d) /  compiled  by  Maxwell  in 
the  late  19th  century.  They  describe  the  interrelation  of 
the  charge  density,  p,  current,  J,  electric  field,  E,  and 
magnetic  field,  B,  (Griffiths  1989) : 


Gauss'  law 

V*E  =  -^ 

(3.1a) 

^0 

(no  name) 

V»B  =  0 

(3.1b) 

Faraday's  law 

VxE  =  -:^ 
dt 

(3.1c) 

Ampere ' s  law 

VxB  =  ;U„J  +  //o£o  — 
at 

(3. Id) 

Maxwell's  equations  provide  a  conceptual  picture  of  the 
fields  we  should  expect  from  a  return  stroke.  If  we 
consider  a  current  flowing  in  a  long  straight  conductor. 
Gauss'  law  predicts  a  change  in  the  electric  field  along  the 
conductor  due  to  the  charge  carried  from  one  end  to  the 
other  by  the  current.  Ampere's  law  predicts  that  the 
current  (and  electric  field  change)  will  produce  a  magnetic 
field  in  the  azimuthal  direction  (in  accordance  with  the 
right-hand  rule) .  If  the  conductor  is  vertical  (as  we 
approximate  a  return  stroke  channel  to  be) ,  the  electric 
field  change  will  also  be  vertical,  and  the  magnetic  field 
will  be  azimuthal.  This  is  why  the  magnetic  direction 
finders  discussed  in  Chapter  II  sense  the  vertical  electric 
field  component  and  the  horizontal  magnetic  field 
components . 

b.  Solution  to  Maxwell 's  equations 

The  first  step  in  obtaining  a  solution  to  Maxwell's 
equations  is  to  formulate  the  scalar  and  vector  potentials, 

(j)  and  A  (Griffiths  1989;  Oman  1987;  Uman  et  al.  1975)  ; 
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The  potentials  are  related  by  (Griffiths  1989) : 

di 

V*K+fioeo  —  =  0  (3.4) 

The  position  vectors  r,  r',  and  R/  are  related  by: 

R  =  r-r'  (3.5) 

The  electric  and  magnetic  fields  expressed  in  terms  of 
the  scalar  and  vector  potentials  are  (Griffiths  1989) : 

E(r,0  =  -V<^-—  (3.6) 

B(r,0  =  VxA  (3.7) 

We  solve  by  substituting  for  the  potentials  in  these 
expressions.  The  boundary  condition  on  a  conductor  dictates 
that  there  be  no  electric  field  component  tangent  to  the 
conductor.  Thus,  if  we  think  of  the  ground  as  a  perfect 
conductor,  the  electric  field  must  have  only  a  vertical 
component,  Ez,  and  the  magnetic  field  must  have  only  an 
azimuthal  component,  B^.  The  solution  in  this  case  is 
obtained  using  the  "method  of  images",  whereby  the  fields 
are  computed  including  mirror  images  charges  and  currents  of 
opposite  sign  on  the  opposite  side  of  the  conductor.  This 
produces  an  electric  field  with  only  a  component  normal  to 
the  surface. 

Using  the  method  of  images  to  compute  the  potentials 
and  substituting  into  equations  3.6  and  3.7,  we  obtain  the 
following  solutions  for  the  vertical  electric  field  and  the 
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azimuthal  magnetic  flux  density  (based  on  Uman  1987  and  Uman 
et  al.  1975): 


EXrJ) 
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(3.8) 


Bt(r,t)  = 


/fo. 

Iti 


\  r  , 

't 

\  r  a(h, t-RIc) 


f 


dh 


(3.9) 


cR^  a 

If  we  assume  that  r  is  large  compared  to  h  (i.e.,  we  are 
a  long  distance  from  the  channel),  we  can  neglect  h  and  make 
the  following  approximation: 
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cr  a 

The  three  electric  field  terms  are:  the  electrostatic  term, 
proportional  to  r~^ ,  produced  by  the  time  integral  of 
current;  the  induction  term,  proportional  to  r~^ ,  produced  by 


the  channel  current;  and  the  radiation  term,  proportional  to 
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r~^ ,  produced  by  the  time  derivative  of  the  current.  The  two 
magnetic  field  terms  are  the  induction  and  radiation  terms, 
analogous  to  their  electric  field  counterparts  (Uman  1987) . 

Near  the  radiation  source  all  three  electric  field 
terms  and  both  magnetic  field  terms  contribute  significantly 
to  the  total  radiation,  complicating  the  solution 
considerably.  But  since  r~^  and  r~^  approach  0  much  more 
rapidly  than  r~^ ,  we  can  neglect  the  electrostatic, 
induction,  and  magnetostatic  terms  at  large  distances, 
leaving  only  the  radiation  terms  (Equations  3.12  and  3.13). 
We  can  make  this  so-called  radiation  field,  or  far  field, 
approximation  without  much  loss  of  accuracy  when  r  is  many 
times  the  height  of  the  channel  (Uman  et  al.  -1975;  Uman 
1987;  Griffiths  1989)  . 


EXr,t)^- 


1  t\\  ci{h,t - r ! c) 


Ins^c  r 


i 


dh 


//o  cM  a(h,t-r  ! c) 


dh 


(3.12) 


(3.13) 


Itvc  r  <3 

If  we  make  the  further  assumptions  that  the  return 
stroke's  velocity  up  the  channel,  v,  and  current  waveform 
are  both  constant  in  time  (as  in  the  transmission  line  model 
discussed  in  Chapter  V)  we  can  eliminate  the  integral  and 


the  current  derivative  as  follows: 


d1  dh 

1 = -'’t  Jif = ■'’i! 

The  solutions  then  take  the  form: 

EXr,  t)  «  ^[KKJ  -ric)-  i(h„  i-r/c)]  (3.14) 

2nr 

B^{r, t)  »  -rlc)-i{\,t-rlc)]  (3.15) 

2wcr 

Thus,  the  radiation  field  strength  is  proportional  to  the 
return  stroke  velocity,  v,  and  inversely  proportional  to  the 
distance,  r,  from  the  channel.  The  t-r/c  appearing  in  the 
current  expression  is  the  retarded  time  and  accounts  for  the 
time  required  for  the  radiation  to  travel  from  the  channel 
to  the  observer. 

The  temporal  waveform  of  the  fields  is  initially 
identical  to  the  waveform  at  the  channel  base,  hb.  As  the 
return  stroke  wavefront  reaches  the  top  of  the  channel,  the 
field  waveform  becomes  a  "mirror  image"  of  the  current  at 
the  channel  top,  /»,.  These  "turn-on"  and  "turn-off"  fields 
have  been  documented  in  the  literature  (Uman  et  al .  1975  and 
LeVine  and  Meneghini  1978)  but  generally  are  not  observed  in 
nature,  since  we  do  not  expect  an  actual  return  stroke 
channel  to  behave  as  a  straight  vertical  segment. 
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LeVine  and  Meneghini  (1978)  point  out,  however,  that 
the  turn-on  and  turn-off  effects  are  the  source  of  the 
noisiness  of  actual  return  stroke  radiation.  They  suggest 
that  the  channel  can  be  represented  by  a  series  of  randomly 
oriented  segments  connected  at  their  endpoints .  The  channel 
then  radiates  from  each  junction  point  between  segments 
where  the  return  stroke  current  turns  on  and  off.  The 
intensity  of  the  radiation  detected  from  each  junction  point 
depends  on  the  adjoining  segments'  orientation  relative  to 
each  other  and  to  the  observer  (LeVine  and  Meneghini  1978) . 
b.  Documented  return  stroke  electric  and  magnetic  fields 

Return  stroke  radiation  has  been  analyzed  extensively 
in  the  literature  using  a  variety  of  sensors  with  a  range  of 
capabilities  (Uman  et  al.  1970;  Weidman  and  Krider  1978;  Lin 
et  al.  1979;  Krider  et  al.  1980;  Willet  et  al.  1989;  Cooray 
and  Perez  1994) .  Typical  characteristics  of  documented 
return  stroke  fields  are  summarized  here.  Examples  of 
measured  return  stroke  radiation  are  depicted  in  Fig.  3.1. 

Signal  amplitudes  are  typically  normalized  to  a  range 
of  100  km  in  the  literature;  the  NLDN  uses  this  convention 
as  well.  This  is  accomplished  (assuming  the  r~^  signal 
strength  range  dependence)  by  multiplying  the  measured 
signal  strength  by  a  factor  of  (r/ 100  km)  (Orville  1991; 

Idone  et  al.  1993;  Hiscox  and  Cummins  1993) .  Typical 
normalized  return  stroke  peak  electric  fields  are  6  to  8  V 
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Fig.  3.1.  Typical  first  (bottom)  and  subsequent  return 
stroke  waveforms  with  dart-stepped  leaders  (middle)  and  dart 
leaders  (top)  (adapted  from  Uman  1987) . 
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for  first  return  strokes  and  4  to  6  V  m“^  for  subsequent 
return  strokes  (Uman  1987) . 

Return  stroke  fields  typically  rise  to  peak  strength  in 
two  stages:  an  initial  "slow  front"  lasting  a  few  ^.s  (Uman 
1987;  Weidman  and  Krider  1978);  and  a  "fast  transition"  in 
which  the  signal  rises  from  about  10%  to  90%  of  peak  signal 
strength  in  approximately  .1  ^s  (Uman  1987) .  Subsequent 
strokes  have  faster  slow  front  risetimes  (.5-1  |is)  and  do 
not  reach  peaks  as  high  as  first  strokes  (Weidman  and  Krider 
1978) . 

Pulsewidth  after  the  initial  peak  is  typically  tens  of 
US  (Uman  1987)  .  After  the  peak,  signal  strength  drops  off 
more  slowly  and  the  waveform  may  have  a  "shoulder"  of  a  few 
US  width.  First  return  strokes  have  multiple  subsidiary 
peaks  at  10  to  30  (j,s  intervals  corresponding  to  current 
flowing  in  the  channel  branches  (Weidman  and  Krider  1978) . 

As  discussed  in  Chapter  II,  the  wide  range  of 
timescales  in  lightning  flash  processes  causes  radiation 
over  a  wide  frequency  range.  Higher  frequencies  are 
attenuated  preferentially  in  propagation  over  a  surface  with 
finite  conductivity,  which  creates  uncertainty  in  the  actual 
upper  frequency  limit  of  return  stroke  radiation.  However, 
documented  return  stroke  radiation  has  been  measured  well 
into  the  hundreds  of  MHz  (Uman  1987) . 
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2.  Peak  Current  Calibration 

Since  the  radiation  approximation  discussed  above 
results  in  field  strengths  that  are  proportional  to  the 
channel  current,  we  can  solve  the  reverse  problem — 
determining  channel  current  from  measured  fields — by 
rearranging  the  expression: 

2/r;  In 

W  «  - ^  (3.16) 

This  linear  relation  between  peak  field  strength  and  peak 
current  is  the  basis  for  estimating  return  stroke  peak 
current  from  remote  field  measurements  (Rachidi  and 
Thottappillil  1992;  Rakov  et  al.  1992) .  Note  that  this 
requires  knowledge  of  the  return  stroke  velocity,  v,  which 
we  have  assumed  to  be  constant.  Since  there  is  no  way  of 
determining  v  for  each  detected  stroke  we  cannot  obtain  a 
theoretical  conversion  constant  from  field  strength  to  peak 
current;  a  calibration  must  be  accomplished  to  determine  it 
empirically. 

a.  Orville  calibration 

Orville  (1991)  was  the  first  to  publish  such  a 
calibration  for  a  network  of  magnetic  DFs.  He  accomplished 
it  using  18  return  strokes  triggered  by  a  French  research 
group  at  the  Kennedy  Space  Center  during  1985-88.  The 
current  of  the  triggered  return  strokes  was  measured 
directly  at  the  channel  base  and  their  radiation  was 
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detected  by  6  DFs  from  the  SUNYA  network  in  Florida  and 
Georgia  with  known  distances  from  the  triggering  site. 

He  computed  the  average  range-normalized  signal 
strength  {RNSS)  from  each  DF  that  detected  a  triggered 
return  stroke,  then  derived  the  following  linear  best-fit 
relation  between  the  average  RNSS  and  the  measured  peak 
current : 

I==2.3kA  +  (.I9kALLF')RNSS  (3.17) 

The  standard  deviation  of  the  data  points  about  the  best-fit 
line  is  approximately  6  kA  (Fig.  3.1). 

This  work  validated  the  approximate  rule  of  thumb  used 
to  date  with  LLP  DFs:  dividing  RNSS  by  5  to  obtain  peak 
current  (Hiscox  and  Cummins  1993) .  This  rule  of  thumb 
produces  results  within  approximately  6%  of  Orville's 
calibration. 

As  a  further  step,  Orville  calibrated  for  each  DF 
individually  and  determined  that  measured  signal  strength 
fell  off  as  due  to  attenuation. 

b.  Idone  et  al.  calibration 

In  1993,  Idone  et  al.  published  a  reexamination  of 
Orville's  calibration,  using  an  expanded  data  set  of  57 
triggered  return  strokes  from  36  different  flashes.  Their 
best-fit  relation,  using  the  same  technique  as  Orville,  is: 


I  =  4.20kA  +  (171  kA  LLF')RNSS 


(3.18) 
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with  a  comparable  standard  deviation. 

If  Orville's  power  law  relationship  is  used  to 

compute  the  normalized  signal  strength  {RNSS=(r/J00  ,  the 

best-fit  relation  becomes: 

I^5.20kA  + (.145kALLF^)RNSS  (3.19) 

This  is  the  calibration  adopted  by  Global  Atmospherics  for 
use  with  the  NLDN  (Hiscox  and  Cummins  1993) .  Idone  et  al. 
estimate  a  peak  current  uncertainty  of  10-15%  in  flashes 
with  peak  currents  between  15  and  60  kA  (Fig.  3.2). 
c.  Forcing  calibration  best-fit  intercept  to  0  kA 

The  non-zero  intercepts  in  the  Orville  and  Idone  et  al. 
calibrations  are  troubling  and  have  yet  to  be  explained. 
Based  on  Eq.  3.16,  we  expect  the  best-fit  line  relating  peak 
current  and  RNSS  to  pass  through  the  origin. 

As  a  result  of  the  increased  number  of  weak  flashes 
detected  by  the  NLDN  since  the  upgrade  (discussed  in  Chapter 
IV) ,  Global  Atmospherics  has  been  forced  to  deal  with  this 
issue.  It  has  introduced  a  new  peak  current  calibration 
with  a  0  kA  intercept,  I  =  (.185  kA  LLF‘)RNSS  ,  effective  May  1996 
(C\ammins  et  al.  1996)  . 

Best-fit  relations  using  the  same  data  as  Orville  and 
Idone  et  al.  but  with  intercepts  constrained  to  0  kA  results 
in  higher  slopes  (resulting  in  larger  estimated  peak 
currents)  and  increased  standard  deviations  of  the  data 
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about  the  lines  (since  the  relation  is  no  longer  a  best-fit 
to  the  data) .  The  slope  of  the  Orville  calibration 
increases  from  .178  to  .195  kA  LLP“^  and  the  standard 
deviation  increases  from  6.304  to  6.427  kA;  the  Idone  et  al. 
calibration  slope  increases  from  .171  to  .203  kA  LLP"^  and 
the  standard  deviation  increases  from  4.649  to  4.984  kA. 

Appendix  A  contains  a  verification  of  the  Orville  and 
Idone  et  al.  calibrations  using  their  published  raw  data,  a 
correction  to  Orville's  calibration  to  account  for  an  error 
in  one  DF  location,  and  the  recalibration  of  both  data  sets 
with  forced  0  kA  intercepts. 
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CHAPTER  IV 
DATA  ANALYSIS 

The  purpose  of  this  chapter  is  to  describe 
quantitatively  the  change  in  mean  peak  current,  flash  count, 
and  percent  positive  flashes  observed  since  the  NLDN 
upgrade . 

1 .  Previous  work 

In  a  study  of  thunderstorm  prediction  in  the  Whiteman 
Air  Force  Base,  Missouri  vicinity  Bass  (1996)  documented 
several  changes  in  NLDN  archived  data:  1)  the  percent  of 
flashes  in  1994  and  1995  having  positive  polarity  doiibled 
compared  to  1989-93,  2)  the  1995  mean  peak  current  decreased 
10-15%  for  negative  flashes  and  20%  for  positive  flashes, 
and  3)  in  1995  twice  as  many  negative  flashes  with  peak 
currents  less  than  20  kA  were  detected  and  most  of  the 
positive  flash  count  increase  was  for  flashes  with  peak 
currents  less  than  30  kA. 

2.  Data  set  and  analysis  methods 

a.  1989-95  NLDN  archived  data 

The  NLDN  records  the  time,  location,  polarity,  RNSS, 
and  multiplicity  of  every  CG  lightning  flash  it  detects. 
These  data  are  encoded  in  binary  form,  using  eleven  bytes 
per  flash  and  disseminated  to  NLDN  subscribers  in  near  real¬ 
time.  This  analysis  uses  NLDN  lightning  data  archived  from 
1989  through  1995. 
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Flash  counts  and  mean  peak  currents  are  analyzed  on  a 
monthly  and  yearly  basis  for  both  positive  and  negative 
flashes.  To  provide  for  regional  comparisons  and  to  develop 
the  initial  data  analysis  process,  three  regions  within  the 
US  data  set  were  chosen  to  analyze  before  undertaking  the 
whole  archive.  Table  4.1  lists  the  three  regions  and  their 
boundaries.  Each  contains  approximately  ten  million  flashes 
for  the  seven  year  period. 


Table  4.1.  Data  analysis  geographical  subsets 


Area 

Latitude  bounds  (°N) 

Longitude  bounds (°W) 

Kansas 

36.0  -  41.0 

93.25  -  103.5 

Ohio  Valley 

36.0  -  41.0 

78.0  -  88.25 

Florida 

27.0  -  29.5 

81.25  -  82.25 

After  completing  the  analysis  for  each  region,  the  full 
US  data  set  of  over  134  million  flashes  was  analyzed  in 
greater  detail.  Positive  and  negative  monthly  and  yearly 
flash  counts  and  mean  peak  currents  were  computed  as  in  the 
original  three  regions.  In  addition,  flash  count 
distributions  as  a  function  of  peak  current  were  computed  to 
determine  if  the  increase  in  flash  count  occurred 
preferentially  for  stronger  or  weaker  flashes. 
b.  Analysis  tools 

The  NLDN  archived  data  were  initially  analyzed  with  the 
FLASH  lightning  analysis  software,  a  DOS-based  package  that 
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uses  a  monochrome  monitor  as  a  user  interface  and  a  color 
monitor  to  display  graphical  results.  FLASH  can  be  used  to 
produce  scatter  plots,  time  series,  and  contours  of  flash 
counts,  peak  currents,  multiplicity,  polarity,  and  other 
parameters  of  detected  CG  flashes.  It  reads  NLDN  data  in 
the  binary  format  described  above. 

FLASH  produces  useful  results,  but  can  be  very  time 
consuming  to  use  because  it  has  no  capability  for  running 
scripts,  so  the  user  must  input  commands  and  record  results 
manually.  For  this  reason,  the  author  wrote  FORTRAN  code 
that  reads  binary  NLDN  data  files,  computes  monthly, 
seasonal,  and  yearly  flash  counts,  mean  peak  currents  and 
peak  current  distributions,  and  stores  the  results  in  text 
files . 

The  text  files  containing  flash  counts,  mean  peak 
currents,  and  peak  current  distributions  are  then  read  with 
the  MathCAD  6.0+  software  package,  where  further  statistical 
analysis  was  performed  (i.e.  computed  weighted  mean  peak 
currents  and  standard  deviations)  and  preliminary  graphical 
results  were  produced.  Final  graphical  results  were 
produced  using  Axum  5.0  software.  A  sample  of  the 
statistical  analysis  is  contained  in  Appendix  B. 

3.  Results 

Tables  4.2  and  4.3  summarize  the  flash  counts  and  mean 
peak  currents  for  each  of  the  subset  regions  and  for  the 


38 


Table  4.2.  Post-  vs.  pre-upgrade  negative  flash  count  and 
mean  peak  current  comparison 


Kansas 

Ohio 

Florida 

US 

Flash  count 

1989-93  mean 

flash  count  (xlO®) 

1.199 

1.409 

1.311 

16.65 

1994  flash  count  (xlO®) 

Z 

1.692 

1.16 

1.869 

1.03 

1.531 

1.23 

23.02 

1.87 

1995  flash  count  (xlO®) 

Z 

1.144 

-0.13 

1.796 

0.86 

1.195 

-0.65 

20.61 

1.16 

Mean  peak 

current 

1989-93  mean  peak 
current  (kA) 

32.9 

33.7 

41.2 

37.5 

1994  mean  peak  current 

Z 

29.8 

-3.36 

31.8 

-1.37 

36.2 

-0.68 

34.26 

1.73 

1995  mean  peak  current 

Z 

26.1 

-7.26 

26.5 

-5.07 

33.2 

-1.08 

30.2 

-3.39 

Z  =  standard  deviations  above  pre-upgrade  mean 
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Table  4.3.  Post-  vs.  pre-upgrade  positive  flash  count  and 
mean  peak  current  comparison 


Kansas 

Ohio 

Florida 

US 

Flash  count 

1989-93  mean 

flash  count  (xlO®) 

0.079 

0.037 

0.038 

0.696 

1994  flash  count  (xlO®) 

Z 

0.138 

1.86 

0.052 

1.02 

0.061 

1.85 

1.193 

2.19 

1995  flash  count  (xlO®) 

Z 

0.148 

2.20 

0.164 

8.57 

0.128 

7.24 

2.107 

6.22 

Mean  peak 

current 

1989-93  mean  peak 
current  (kA) 

55.2 

54.2 

33.9 

54.4 

1994  mean  peak  current 

Z 

44.7 

-8.64 

48.3 

-0.72 

26.5 

-1.11 

45.8 

-1.87 

1995  mean  peak  current 

Z 

37.6 

-14.53 

23.4 

-3.79 

19.3 

-2.21 

31.6 

-4.97 

Z  =  standard  deviations  above  pre-upgrade  mean 
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entire  US  before  the  upgrade  (1989-93),  during  the  upgrade 
(1994),  and  after  the  upgrade  (1995).  For  purposes  of 
comparison,  1989-93  are  considered  pre-upgrade  years  and 
1995  is  considered  the  post-upgrade  year.  Since  the  network 
upgrade  occurred  during  1994,  it  is  considered  separately 
from  the  pre-  and  post-upgrade  data  sets, 
a.  Kansas 

The  Kansas  data  set  consisted  of  slightly  over  9.5 
million  flashes,  approximately  8.8  million  negative  and  680 
000  positive: 

1)  Negative  flashes 

As  Table  4.2  and  Fig.  4.1  show,  the  post-upgrade 
negative  flash  count  remains  close  to  the  pre-upgrade  mean 
of  approximately  1.2  million  flashes  per  year.  1994  and 
1995  monthly  flash  counts  are,  with  three  exceptions,  within 
two  standard  deviations  of  the  pre-upgrade  means.  The  1995 
negative  flash  count  is  approximately  1.14  million  flashes, 
only  .13  standard  deviation  below  the  1989-93  Kansas  mean 
negative  flash  count. 

While  the  negative  flash  count  remains  relatively 
constant,  the  Kansas  negative  mean  peak  current  decreases 
significantly  after  the  upgrade,  as  depicted  in  Table  4.2 
and  Fig.  4.2.  From  a  pre-upgrade  value  of  32.9  kA  the  1994 
mean  peak  current  drops  3.36  standard  deviations  to  29.8  kA 
and  the  1995  mean  sinks  to  26.1  kA,  a  decrease  of  7.26 
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Kansas  Annual  Negative  Rash  Count 
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Fig.  4.1.  Kansas  annual  (top)  and  monthly  (bottom)  negative 
flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 


Mean  Peak  Current  (kA)  Mean  Peak  Current  (kA) 


42 


Fig.  4.2.  Kansas  annual  (top)  and  monthly  (bottom)  mean 
peak  current,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Kansas  Annual  Positive  Rash  Count 
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Fig.  4.3.  Kansas  annual  (top)  and  monthly  (bottom)  positive 
flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.4.  Kansas  annual  (top)  and  monthly  (bottom)  positive 
mean  peak  current,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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standard  deviations.  Mean  peak  currents  first  exceed  two 
standard  deviations  below  the  pre-upgrade  mean  in  July  1994; 
every  month  of  1995  has  a  mean  peak  current  less  than  two 
standard  deviations  below  the  pre-upgrade  monthly  mean. 

2)  Positive  flashes 

Unlike  the  Kansas  negative  flash  count,  the  positive 
flash  count  increases  significantly  after  the  upgrade,  as 
depicted  in  Table  4.3  and  Fig.  4.3.  From  a  pre-upgrade  mean 
of  approximately  79  000  flashes  per  year,  the  1995  total 
increases  by  2.2  standard  deviations  to  148  000  flashes. 

The  monthly  flash  count  shows  that  this  increase  occurs 
primarily  in  the  late  spring  and  early  summer  of  1995. 

The  Kansas  positive  mean  peak  current  decreases  even 
more  dramatically  than  the  negative  mean  peak  current,  from 
a  pre-upgrade  mean  of  55.2  kA  to  a  1994  mean  of  44.7  kA,  an 
8.64  standard  deviation  decrease,  and  a  1995  mean  of  37.6 
kA,  a  14.53  standard  deviation  decrease.  Fig.  4.4  shows 
that  this  decrease  begins  in  July  1994,  as  it  did  for 
negative  flashes, 
b.  Ohio  Valley 

The  Ohio  Valley  data  set  consists  of  approximately  11.1 
million  flashes,  10.7  million  negative  and  just  under  400 
000  positive: 
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1)  Negative  flashes 

Table  4.2  and  Fig.  4.5  show  that  the  Ohio  Valley 
negative  flash  count  averages  approximately  1.4  million 
flashes  yr"^  before  the  upgrade.  1994  and  1995  both  have 
approximately  1.8  million  flashes,  an  increase  of  less  than 
one  standard  deviation  after  the  upgrade.  As  in  the  Kansas 
data  set,  with  few  exceptions  the  post-upgrade  monthly 
negative  flash  counts  are  all  within  two  standard  deviations 
of  the  pre-upgrade  mean  flash  count.  The  large  flash  count 
in  1993  corresponds  to  the  Midwest  flood  during  that  summer. 
If  1993  were  not  considered,  the  negative  flash  count 
standard  deviation  would  be  much  lower  and  the  post-upgrade 
flash  count  increase  would  be  statistically  more 
significant . 

Ohio  Valley  negative  mean  peak  currents  exhibit  a  very 
similar  pre-  to  post-upgrade  change  as  Kansas  mean  peak 
currents,  as  shown  in  Fig.  4.6.  From  a  pre-upgrade  mean 
peak  current  of  33.7  kA,  the  1994  mean  decreases  to  31.8  kA, 
a  1.37  standard  deviation  increase,  and  the  1995  mean 
decreases  further  to  26.5  kA,  a  5.07  standard  deviation 
decrease.  For  the  Ohio  valley,  the  decrease  in  negative 
mean  peak  current  appears  to  begin  in  September  or  October 


1994. 
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Fig.  4.5.  Ohio  Valley  annual  (top)  and  monthly  (bottom) 
negative  flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.6.  Ohio  Valley  annual  (top)  and  monthly  (bottom) 
negative  mean  peak  current,  1989-95.  Error  bars  are  two 
standard  deviations  above  and  below  the  mean. 


Fig.  4.7.  Ohio  Valley  annual  (top)  and  monthly  (bottom) 
positive  flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.8.  Ohio  Valley  annual  (top)  and  monthly  (bottom) 
positive  mean  peak  current,  1989-95.  Error  bars  are  two 
standard  deviations  above  and  below  the  mean. 


■ 
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2)  Positive  flashes 

As  shown  in  Table  4.3  and  Fig  4.7,  the  Ohio  Valley 
post-upgrade  positive  flash  count  increases  much  more 
dramatically  than  the  Kansas  positive  flash  count — from  a 
pre-upgrade  mean  of  approximately  37  000  flashes  yr"^  to  a 
1994  count  of  52  000,  about  a  one  standard  deviation 
increase,  and  a  1995  count  of  164  000,  an  8.57  standard 
deviation  increase.  As  with  Kansas,  this  appears  to  be 
primarily  a  warm-season  effect,  with  the  1995  increase 
occurring  mainly  during  the  months  April  through  August. 

While  the  Ohio  Valley  positive  mean  peak  current 
decrease  is  more  modest  than  Kansas'  in  terms  of  standard 
deviations  below  the  pre-upgrade  mean,  the  actual  numerical 
change  is  larger.  This  is  due  to  the  large  variance  in  the 
Ohio  Valley  data  caused  by  the  abnormally  high  positive  mean 
peak  current  of  70.6  kA  in  1989.  The  pre-upgrade  mean  peak 
current  is  54.2  kA;  in  1994  the  mean  drops  less  than  one 
standard  deviation  (but  6  kA)  to  48,3  kA  and  in  1995  it 
decreases  further  to  23.4  kA,  only  3.79  standard  deviations, 
but  over  50%  less  than  the  pre-upgrade  value.  Fig.  4.8 
shows  the  decrease  to  have  begun  in  about  November  1994. 
c .  FI ori da 

The  Florida  data  set  consists  of  approximately  9.7 
million  flashes,  9.3  million  negative  and  just  under  380  000 
positive : 
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1)  Negative  flashes 

As  depicted  in  Table  4.2  and  Fig.  4.9,  Florida  negative 
flash  counts  average  just  over  1.3  million  flashes  yr”^ 
before  the  upgrade.  In  1994  the  count  increases  1.23 
standard  deviations  to  over  1.5  million,  but  in  1995  it 
decreases  to  just  under  1.2  million.  Every  1994  and  1995 
monthly  negative  flash  count  except  one  is  within  two 
standard  deviations  of  the  pre-upgrade  monthly  negative  mean 
flash  count. 

As  is  the  case  with  positive  flashes  in  the  Ohio 
Valley,  one  abnormally  high  annual  peak  current  (1990) 
causes  a  large  variance  in  the  pre-upgrade  Florida  negative 
mean  peak  current,  so  a  relatively  large  post-upgrade  change 
produces  less  meaningful  statistical  results.  The  pre¬ 
upgrade  negative  mean  peak  current  is  41.2  kA,  approximately 
8  kA  higher  than  the  Kansas  and  Ohio  Valley  negative  mean 
peak  currents.  In  1994  this  decreases  to  36.2  kA,  a  12% 
decrease,  but  only  .68  standard  deviations;  in  1995  the  mean 
peak  current  sinks  to  33.2  kA,  a  19%  decrease  from  the  pre¬ 
upgrade  mean,  but  only  1.08  standard  deviations.  As  Fig. 

4.10  shows,  monthly  mean  peak  currents  after  April  1994  are 
below  the  pre-upgrade  mean,  but  relatively  few  exceed  two 


standard  deviations  below  the  mean. 
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Fig.  4.9.  Florida  annual  (top)  and  monthly  (bottom) 
negative  flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.10.  Florida  annual  (top)  and  monthly  (bottom) 
negative  mean  peak  current,  1989-95.  Error  bars  are  two 
standard  deviations  above  and  below  the  mean. 
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Fig.  4.11.  Florida  annual  (top)  and  monthly  (bottom) 
positive  flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.12.  Florida  annual  (top)  and  monthly  (bottom) 
positive  mean  pealc  current,  1989-95.  Error  bars  are  two 
standard  deviations  above  and  below  the  mean. 
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2)  Positive  flashes 

As  in  the  Kansas  and  Ohio  Valley  data,  the  Florida 
positive  flash  count  increases  more  dramatically  after  the 
upgrade  than  the  negative  flash  count.  There  are  an  average 
of  just  over  38  000  flashes  yr"^  before  the  upgrade.  In 
1994,  the  count  increases  by  1.85  standard  deviations  to  61 
000  flashes;  in  1995  it  increases  to  128  000  flashes — a  7.24 
standard  deviation  increase.  Fig.  4.11  shows  that,  as  with 
Ohio  positive  flashes,  the  increase  is  most  dramatic  during 
the  warm  season,  from  May  to  September  1995. 

The  large  variance  that  affected  Florida  negative  mean 
peak  currents  also  affects  Florida  positive  mean  peak 
currents,  causing  large  numerical  changes  to  have  relatively 
little  statistical  significance.  The  pre-upgrade  positive 
mean  peak  current  is  33.9  kA,  about  20  kA  lower  than  the 
Kansas  and  Ohio  Valley  means.  In  1994  the  mean  drops  to 
26.5  kA,  a  1.11  standard  deviation  decrease,  and  in  1995  it 
sinks  to  19.3  kA,  a  2.21  standard  deviation  decrease 
(although  the  value  is  43%  lower) .  The  decrease  is  evident 
in  Fig.  4.12  after  September  1994. 
d.  Continental  United  States 

The  full  US  NLDN  archived  data  set  contains  roughly 
133.7  million  flashes,  126.9  million  of  which  are  negative 
and  6.8  million  of  which  are  positive.  Many  of  the  same 
pre-  and  post-upgrade  flash  characteristics  are  exhibited  by 
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the  full  data  set  as  in  the  smaller  regions.  Due  to  its 
much  larger  sample  size,  though,  it  tends  to  have  less 
"noisiness"  than  the  smaller  samples  in  the  subsets, 
resulting  in  smaller  standard  deviations  for  flash  counts 
and  mean  peak  currents. 

1)  Negative  flashes 

As  Table  4.2  and  Fig.  4.13  show,  the  pre-upgrade  US 
negative  mean  flash  count  is  16.65  million  flashes  yr"\  The 
1994  flash  count  is  1.87  standard  deviations  higher  at  23.02 
million  flashes,  while  the  1995  flash  count  is  only  1.16 
standard  deviations  higher  at  20.61  million  flashes.  As 
with  the  Kansas  and  Ohio  Valley  data  sets,  only  a  handful  of 
post-upgrade  monthly  negative  flash  counts  exceeds  two 
standard  deviations  above  the  pre-upgrade  mean  negative 
flash  count. 

Table  4.2  and  Fig.  4.14  depict  US  negative  mean  peak 
currents.  The  pre-upgrade  mean  is  37.5  kA,  roughly  4  kA 
higher  than  the  Kansas  and  Ohio  Valley  pre-upgrade  means  and 
4  kA  below  the  Florida  pre-upgrade  mean.  In  1994  the  mean 
decreases  1.73  standard  deviations  to  34.3  kA  and  in  1995  it 
decreased  to  30.2  kA,  3.39  standard  deviations  below  the 
pre-upgrade  mean.  The  monthly  negative  mean  peak  current 
first  exceeds  two  standard  deviations  below  the  pre-upgrade 
mean  after  September  1994  and  remains  well  below  that  value 


for  all  of  1995. 
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Fig.  4.13.  US  annual  (top)  and  monthly  (bottom)  negative 
flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.14.  US  annual  (top)  and  monthly  (bottom)  negative 
mean  peaJc  current,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.15.  US  annual  (top)  and  monthly  (bottom)  positive 
flash  count,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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Fig.  4.16.  US  annual  (top)  and  monthly  (bottom)  positive 
mean  pealc  current,  1989-95.  Error  bars  are  two  standard 
deviations  above  and  below  the  mean. 
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2)  Positive  flashes 

While  the  US  negative  flash  count  increases  only 
slightly,  the  positive  flash  count  increases  significantly, 
as  it  does  in  each  of  the  three  regions.  Table  4.3  and  Fig. 
4.15  show  that  the  pre-upgrade  mean  positive  flash  count  is 
approximately  696  000  flashes  yr”^.  In  1994  this  increases 
2.19  standard  deviations  to  1.193  million  flashes.  The  1995 
positive  flash  count  is  6.22  standard  deviations  higher  than 
the  pre-upgrade  mean  at  2.107  million  flashes,  a  threefold 
increase.  Every  month's  positive  flash  count  exceeds  two 
standard  deviations  above  the  pre-upgrade  mean  after 
September  1994. 

Table  4.3  and  Fig.  4.16  depict  US  positive  mean  peak 
currents.  The  pre-upgrade  mean  is  54.4  kA,  in  close 
agreement  with  the  Kansas  and  Ohio  Valley  data  sets,  but 
significantly  higher  than  the  Florida  pre-upgrade  mean.  The 
1994  US  positive  mean  peak  current  is  1.87  standard 
deviations  lower  at  45.8  kA  and  the  1995  mean  is  4.97 
standard  deviations  lower  at  31.6  kA.  As  for  negative 
flashes,  September  1994  marks  the  point  where  monthly 
positive  mean  peak  currents  first  exceed  two  standard 
deviations  below  the  pre-upgrade  monthly  mean  peak  currents, 
e.  US  peak  current  distributions 

We  conclude  the  NLDN  data  analysis  by  looking  at  the 
distribution  of  peak  currents  in  negative  and  positive 
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flashes.  The  purpose  for  doing  this  is  to  determine  if  the 
increase  in  NLDN  detected  flash  counts  is  a  dependent  upon 
peak  current. 

1)  Negative  flashes 

Fig.  4.17  is  a  histogram  of  US  annual  negative  flash 
counts  binned  in  1  kA  intervals.  After  the  upgrade,  the 
mode  of  the  negative  peak  current  distribution  decreases 
from  pre-upgrade  values  of  between  24  and  26  kA  to  20  kA  in 
1995.  The  flash  count  also  increases  slightly,  as  discussed 
earlier.  The  bottom  half  of  Fig.  4.17  depicts  the  same  data 
but  with  the  pre-upgrade  years  (1989-93)  replaced  by  the 
pre-upgrade  mean  flash  count  distribution  with  error  bars 
two  standard  deviations  above  and  below. 

Fig.  4.18  depicts  the  same  data  normalized  to  the  1989- 
93  mean  annual  flash  count  distribution.  The  top  half  of 
the  figure  shows  the  ratio  of  flashes  detected  to  the  pre¬ 
upgrade  mean  flash  count  at  each  peak  current  interval.  The 
bottom  half  depicts  the  same  data  in  terms  of  standard 
deviations  above  or  below  the  mean.  This  figure  clearly 
shows  that  the  post-upgrade  increase  in  negative  flash  count 
is  confined  to  flashes  having  peak  currents  less  than  30  kA. 
Above  30  kA  the  1995  flash  counts  are  actually  below  the 
pre-upgrade  mean.  By  contrast,  7  times  the  pre-upgrade  mean 
number  of  flashes  is  detected  at  5  kA  in  1995.  At  6  kA,  the 


1995  flash  count  is  nineteen  standard  deviations  above  the 
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Fig.  4.17.  US  negative  peak  current  distribution,  1989-95. 
Pre-upgrade  years  are  depicted  individually  at  top  and  the 
pre-upgrade  mean  is  depicted  at  bottom.  Error  bars  are  two 
standard  deviations  above  and  below  the  1989-93  mean. 
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Ratio  of  1994-5  Negative  Fiash  Count  to  Pre-Upgrade  Fiash  Count 


1994-95  Negative  Fiash  Count  Standard  Deviations  Above  1989-93  Mean 


Fig.  4.18.  US  negative  pre-  to  post-upgrade  flash  count 
ratio.  1994  and  '95  flash  counts  are  expressed  as  a  ratio 
of  the  pre-upgrade  mean  at  top  and  in  terms  of  standard 
deviations  above  the  pre-upgrade  mean  at  bottom. 
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Fig.  4.19.  US  positive  peak  current  distribution,  1989-95. 
Pre-upgrade  years  are  depicted  individually  at  top  and  the 
pre-upgrade  mean  is  depicted  at  bottom.  Error  bars  are  two 
standard  deviations  above  and  below  the  1989-93  mean. 
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Ratio  of  1994-5  Positive  Flash  Count  to  Pre-Upgrade  Fiash  Count 


1994-95  Positive  Fiash  Count  Standard  Deviations  Above  1989-93  Mean 


Fig.  4.20.  US  positive  pre-  to  post-upgrade  flash  count 
ratio.  1994  and  '95  flash  counts  are  expressed  as  a  ratio 
of  the  pre-upgrade  mean  at  top  and  in  terms  of  standard 
deviations  above  the  pre-upgrade  mean  at  bottom. 
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pre-upgrade  mean  flash  count.  At  all  peak  currents  less 
than  24  kA  the  1995  flash  count  is  in  excess  of  two  standard 
deviations  greater  than  the  pre-upgrade  mean  flash  count. 

2)  Positive  flashes 

Like  the  negative  flash  count  distribution,  the 
positive  flash  count  distribution  is  right-tailed,  as  shown 
in  Fig.  4.19.  Its  mode  is  between  12  and  14  kA  during  1989- 
93;  the  1995  mode  is  also  12  kA.  The  post-upgrade  increase 
in  flash  count  is  more  evident  than  in  the  negative  peak 
current  distribution. 

Fig.  4.20  depicts  the  positive  flash  count  distribution 
normalized  to  the  1989-93  mean.  The  increase  in  weak  flash 
counts  is  even  more  pronounced  than  for  negative  flashes. 

The  1995  2  kA  positive  flash  count  is  58  times  greater  than 
the  pre-upgrade  mean  flash  count,  a  50  standard  deviation 
increase.  At  4  kA,  the  1995  positive  flash  count  is  72 
standard  deviations  higher  than  the  pre-upgrade  mean  flash 
count.  As  With  negative  flashes,  the  increase  is  confined 
to  weaker  flashes;  above  approximately  60  kA,  there  is  no 
appreciable  increase  in  positive  flash  count  after  the 
upgrade,  while  all  flash  counts  below  56  kA  are  increased  by 
more  than  two  standard  deviations. 
f.  Summary 

The  statistical  analysis  of  the  US  NLDN  archived  data 
set  contained  in  this  chapter  has  shown  us  that: 
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1)  The  post-upgrade  negative  flash  count,  while 
increased,  is  not  statistically  different  from  the  pre¬ 
upgrade  value. 

2)  The  post-upgrade  negative  mean  peak  current 
decreases  significantly  after  September  1994. 

3)  Unlike  the  post-upgrade  negative  flash  count,  the 
positive  flash  count  increases  significantly  after  September 
1994. 

4)  The  post-upgrade  positive  mean  peak  current 
decreases  significantly  after  September  1994. 

5)  September  1994  marks  the  point  where  the  upgrade 
first  has  a  statistically  significant  (defined  as  a  change 
in  excess  of  two  standard  deviations)  effect  on  all  observed 
parameters  except  negative  flash  count. 

6)  The  increase  in  both  negative  and  positive  flash 
counts  occurs  preferentially  in  weak  flashes. 
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CHAPTER  V 

RETUEIN  STROKE  DETECTION  MODELING 

The  purpose  of  this  chapter  is  to  model  the  changes 
made  to  the  NLDN  DFs  during  the  1994  network  upgrade  and 
compare  the  resulting  changes  in  flash  counts  and  mean  peak 
currents  to  the  results  from  Chapter  IV.  Two  models  are 
utilized.  The  first  is  a  return  stroke  radiation  model 
which  computes  radiation  waveforms,  peak  signal  strength, 
and  signal  pulsewidth  of  a  model  return  stroke.  The  second 
is  a  detection  model  which  uses  output  from  multiple  runs  of 
the  return  stroke  radiation  model  to  compute  peak  current 
distributions  of  randomly  located  flashes  with  random  peak 
currents  detected  by  DFs.  The  peak  current  distribution  of 
flashes  detected  by  the  model  using  pre-  and  post-upgrade 
pulsewidth  detection  criteria  are  compared  with  pre-  and 
post-upgrade  peak  current  distributions  of  archived  NLDN 
data. 

1 .  Previous  work 

The  processes  composing  a  lightning  flash  have  been 
modeled  in  various  levels  of  detail  for  several  decades. 

Most  of  this  work  has  centered  around  modeling  the 
electromagnetic  radiation  emitted  by  a  return  stroke,  or  the 
opposite  problem,  determining  the  spatial  and  temporal 


structure  of  the  return  stroke  current  from  measured  fields. 
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Over  the  course  of  this  work,  several  return  stroke  current 
models  have  been  developed  and  used, 
a.  Return  stroke  current  models 

All  return  stroke  models  in  the  literature  assume  a 
return  stroke  wavefront  that  travels  up  the  channel,  as  in 
Chapter  II,  with  a  current  waveform  behind  the  wavefront 
that  is  a  function  of  space,  time,  or  both. 

1)  Simple  models 

Uman  et  al.  (1975)  use  a  simple  triangular  wave 
propagating  up  a  channel,  which  has  the  disadvantage  of  two 
discontinuities  in  the  current  derivative.  LeVine  and 
Meneghini  (1978)  avoid  this  difficulty  by  using  a  compound 
exponential  waveform. 

2)  The  Bruce-Golde  model 

The  Bruce-Golde,  or  BG,  model  is  the  oldest  in  use  (it 
dates  back  to  1948)  and  assumes  a  current  that  is  uniform 
behind  the  return  stroke  wavefront  and  0  everywhere  ahead  of 
it  (Uman  and  McLain  1970).  Lin  et  al.  (1980)  showed  that 
this  is  not  a  physically  reasonable  model,  and  the  BG 
model's  use  has  declined  in  recent  years. 

3)  The  transmission  line  model 

The  transmission  line,  or  TL,  and  modified  transmission 
line,  or  MTL,  models  are  a  broad  category  of  models  which 
represent  the  return  stroke  as  a  current  pulse  on  a  lossless 
transmission  line  (Uman  1987) .  This  causes  the  return 


73 


stroke's  waveform  to  remain  constant  in  time  as  it 
propagates  up  the  channel.  TL  models  are  still  widely  used 
in  the  literature  and  will  be  used  in  the  model  discussed  in 
this  chapter. 

4)  The  Lin,  Oman,  and  Standler  model 

In  1980,  Lin  et  al.  introduced  the  most  complex  model 
to  date,  the  LUS  model.  It  has  three  distinct  components: 
a  pulse  component  representing  breakdown  at  the  return 
stroke  wavefront,  a  uniform  component  representing  the 
residual  leader  current,  and  a  corona  component  representing 
the  discharge  of  the  corona  sheath  created  by  the  leader 
process . 

The  LUS  pulse  velocity  is  arbitrary  and  its  risetime  is 
equal  to  that  of  the  measured  fields.  The  uniform  current 
component  is  determined  by  the  measured  vertical  electric 
field  derivative.  The  corona  current  is  modeled  as  a  series 
of  discrete  current  sources  turned  on  as  the  pulse  passes 
and  whose  amplitude  decreases  exponentially  with  height  (Lin 
et  al.  1980) . 

The  corona  and  uniform  current  components  account  for 
most  of  the  charge  transfer  in  the  model.  Channel 
tortuosity,  or  the  degree  of  random  twisting  of  the  channel, 
as  well  as  uncertainty  in  return  stroke  pulse  velocity, 
leads  to  inaccuracies  in  modeled  radiation  (Lin  et  al. 


1980)  . 
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5)  Modifications  to  the  Lin,  Uman,  and  Standler  model 

In  1981,  Master  et  al.  modified  the  LUS  model  by 

introducing  a  pulse  current  amplitude  that  decays  with 
height,  to  reflect  observational  data  (the  MULS  model) .  In 
1994,  Thottappillil  and  Uman  added  a  height-dependent  time 
constant  to  the  channel  discharge  in  the  corona  component 
(the  VDTC  model) . 

6)  Model  evaluations 

Nucci  et  al.  (1990)  evaluate  the  effectiveness  of  all 
of  the  above  models  using  a  common  channel  base  current  for 
each.  Surprisingly,  they  find  all  produced  reasonable 
approximations  to  actual  return  stroke  radiation  fields, 
while  pointing  out  that  no  model  could  adequately  represent 
the  random  noise  imposed  on  return  stroke  radiation  by 
channel  tortuosity  and  branching. 

For  the  specific  purpose  of  computing  peak  fields  from 
peak  currents  and  vice  versa,  Thottappillil  and  Uman  (1993) 
recommend  using  the  TL  model.  During  the  initial  few  (j,s  of 
a  return  stroke,  they  find  it  to  perform  as  well  as  more 
complicated  models.  Since  the  peak  field  typically  occurs 
very  near  the  beginning  of  the  return  stroke  radiation,  the 
TL  model's  loss  of  accuracy  with  time  does  not  affect  peak 
field  calculations. 

For  this  reason,  the  return  stroke  model  used  here  will 
be  a  TL  model  incorporating  a  sharp  initial  peak  traveling 
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up  the  channel  at  constant  velocity  with  an  exponentially 
decaying  current  behind. 
b.  Geometric  effects 

From  Maxwell's  equations,  we  reason  that  the 
orientation  of  the  channel  affects  the  measured  magnitude  of 
the  peak  vertical  electric  field  and  azimuthal  magnetic  flux 
density  (LeVine  and  Meneghini  1978) . 

Due  to  the  "turn-on"  and  "turn-off"  effects  at  each 
channel  segment  junction  discussed  in  Chapter  III,  channel 
tortuosity  results  in  more  temporal  fine  structure  in  the 
return  stroke  radiation  (LeVine  et  al.  1986;  Uman  1987) . 

Channel  branching  profoundly  affects  first  return 
stroke  radiation  fields  by  causing  subsidiary  peaks  as  the 
charge  contained  in  each  branch  flows  to  ground  (Weidman  and 
Krider  1978).  Willet  et  al.  (1995)  measure  18  dB  more 
spectral  power  at  high  frequencies  (from  500  kHz  to  7  MHz) 
in  first  return  strokes  (with  branches)  than  in  subsequent 
return  strokes. 

Most  models,  including  the  one  presented  in  this 
chapter,  consider  neither  tortuosity  nor  branching,  and  so 
the  radiation  they  produce  is  "cleaner"  than  actual  return 
stroke  radiation  signatures,  having  neither  the  high 
frequency  component  nor  the  subsidiary  peaks  characteristic 
of  actual  return  stroke  radiation. 
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c.  Surface  conductivity  effects 

Due  to  the  variability  in  surface  composition,  we 
expect  that  surface  conductivity  should  vary.  In  no  case, 
however,  do  we  expect  the  ground  to  be  a  perfect  conductor. 
Electromagnetic  radiation  propagating  along  a  finitely 
conducting  surface  will  be  attenuated  (Griffiths  1989) . 

This  causes  the  range  normalized  peak  field  strength  of 
distant  signals  to  be  less  than  that  of  close  signals  (Lin 
et  al.  1980)  .  This  is  the  source  of  the  signal 

strength  dependence  found  by  Orville  (1991)  and  Hiscox  and 
Cummins  (1993)  . 

Attenuation  occurs  preferentially  for  high  frequency 
radiation,  causing  the  fine-structure  of  the  return  stroke 
waveform  to  decay  with  distance  from  the  source  (Uman  1987) . 
The  preferential  attenuation  at  high  frequencies  may  also 
contribute  to  a  time  delay  between  components  of  return 
stroke  radiation  having  different  characteristic  frequencies 
(LeVine  et  al.  1986).  Finally,  since  attenuation  is  caused 
by  propagation  along  a  non-perfectly  conducting  surface,  it 
lessens  with  increasing  height  from  which  the  radiation  was 
emitted  (LeVine  et  al.  1986). 

2.  Finite  difference  return  stroke  field  model 

Appendix  C  contains  the  MathCAD  6.0+  file  used  to 
compute  the  radiation  emitted  by  a  model  CG  return  stroke. 
The  user  supplies  channel  geometry,  peak  current,  return 
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stroke  velocity,  and  spatial  and  temporal  waveform.  From 
this  input,  the  model  computes  1)  the  radiated  vertical 
electric  field  strength  and  azimuthal  magnetic  flux  density 
at  a  distance  r  from  the  channel,  2)  the  peak  signal 
strength  and  range  normalized  signal  strength  (in  LLP  units) 
vs.  range,  and  3)  the  signal  pulsewidth  (time  above  the 
detection  signal  strength  threshold)  vs.  range, 
a.  Geometry 

Fig.  5.1  illustrates  the  geometry  of  the  return  stroke 
model.  It  uses  a  straight  vertical  channel  with  specified 
bottom  and  top  heights.  The  observer  is  assumed  to  be  at 
ground  level.  For  the  results  presented  in  this  chapter,  a 
2  km  tall  channel  with  base  at  the  ground  is  used.  The 
channel  has  no  tortuosity  and  no  branches,  so  high  frequency 
components  and  subsidiary  peaks  are  absent  from  the 
radiation. 

The  well-defined  channel  endpoints  result  in  a  temporal 
radiation  field  waveform  that  initially  follows  the  current 
waveform  at  the  bottom  of  the  channel  (due  to  the  "turn  on" 
effect) ,  and  that  becomes  a  mirror  image  of  the  current 
waveform  after  it  reaches  the  top  of  the  channel  (due  to  the 
"turn  off"  effect) .  In  nature,  the  channel  most  often  ends 
in  a  significant  intercloud  segment  with  more  horizontal 
orientation.  This  tends  to  decrease  or  eliminate  the  "turn 


off"  effect,  which  is  not  often  observed  in  nature. 


Fig.  5.1.  Finite  difference  field  model  channel  geometry. 
The  channel  is  2  km  high  and  vertical.  Azimuthal  symmetry 
is  assumed,  so  the  radiation  has  only  a  radial  dependence. 
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b.  Surface  conductivity 

The  field  model  presented  here  treats  the  ground  as  a 
perfect  conductor/  as  do  most  models  in  the  literature. 

This  allows  the  use  of  the  "method  of  images"  discussed  in 
Chapter  III  to  compute  the  fields.  It  also  eliminates  any 
attenuation  effects  on  the  model  results,  so  we  expect  peak 
signal  strength  to  decrease  as  and  the  radiation  fields 
to  retain  their  high  frequency  components. 

c.  Return  stroke  current  model 

The  return  stroke  current  model  presented  here  is  a 
transmission  line  model — the  temporal  waveform  of  the  return 
stroke  current  is  constant.  It  is  a  simplified  version  of 
the  LUS  model,  consisting  of  a  sharp  initial  breakdown  peak 
at  the  return  stroke  wavefront  followed  by  a  current 
decaying  exponentially  with  height  below  it: 

The  wavefront  propagates  up  the  channel  at  a  constant 
velocity,  v,  of  10®  m  s'\  Fig.  5.2  is  a  snapshot  of  the 
channel  current  for  a  30  kA  peak  current  return  stroke  15  |xs 
after  the  return  stroke  initiated  at  the  surface.  Based  on 
conservation  of  charge  in  the  channel,  we  can  use  a 
continuity  equation  to  compute  the  charge  density  as  a 
function  of  height  and  time,  and  find  that  its  shape  is 
identical  to  the  current  waveform.  The  return  stroke 
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I 


Model  Channel  Current,  t=15[xs 


Fig.  5.2.  Model  channel  current  vs.  height  15  |is  after 
return  stroke  initiation.  This  is  a  snapshot  of  the  current 
waveform  for  a  model  return  stroke  with  peak  current  of  30 
kA  (the  peak  falls  between  channel  grid  points  in  this  view, 
causing  the  peak  to  be  truncated  to  27  kA) .  The  wavefront 
propagates  up  the  channel  at  a  velocity  of  10®  m  s"^ . 
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modeled  here  contains  positive  current  flowing  upward — this 
delivers  a  net  positive  charge  upward,  equivalent  to 
delivering  a  net  negative  charge  downward.  Thus,  we  are 
modeling  a  negative  CG  return  stroke. 

d.  Finite  difference  representations  of  field  equations 
We  compute  the  fields  in  the  same  manner  discussed  in 
Chapter  III,  using  the  following  finite  difference 
approximations  for  the  scalar  and  vector  potential 
equations,  making  use  of  the  method  of  images: 
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where  Rj  =  (r^  +  (hj  - zfy  (refer  to  Fig.  5.1). 

The  vertical  electric  field  strength  and  azimuthal 
magnetic  flux  density  are  then  cpmputed  from  the  scalar  and 
vector  potentials: 
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Model  Radiated  Electric  Field  Strength,  i=100km 


CO 


Fig.  5.3.  Electric  field  strength  (top)  and  magnetic  flux 
density  (bottom)  measured  100  km  from  the  30  kA  peak  current 
model  return  stroke  in  Fig.  5.2.  Time  is  measured  from 
initiation  of  the  return  stroke. 
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As  one  would  expect,  the  results  of  the  numerical 
integrations  and  derivatives  are  sensitive  to  changes  in 
grid  sizes,  Ar  and  Ar,  and  time  step  size,  tsi,  all  of  which 
are  adjustable  by  the  user. 

Fig.  5.3  depicts  the  vertical  electric  field  strength, 
Ez,  and  azimuthal  magnetic  flux  density,  measured  a 
distance  of  100  km  from  the  channel  for  the  30  kA  peak 
current  return  stroke  depicted  in  Fig.  5.2.  The  turn-on  and 
turn-off  effects  are  clearly  illustrated.  As  expected  for  a 
negative  CG  flash,  the  initial  vertical  electric  field 
change  is  negative  and  the  initial  azimuthal  magnetic  flux 
density  is  positive  (refer  to  Chapter  III) .  Note  that  Ez  and 
are  roughly  equal  to  each  other  when  expressed  in  LLP 
units.  This  corresponds  to  the  B^  =  Ez/c  proportionality  we 
expect  based  on  the  radiation  field  approximation, 
e.  Signal  strength  and  pulsewidth 

The  model  next  determines  the  peak  signal  strength  vs. 
range  from  the  channel  for  ranges  of  10  km  to  10^  km.  The 
RNSS  vs.  range  is  computed  by  multiplying  the  peak  signal 
strength  by  [r/ 100  km).  The  signal  pulsewidth  vs.  range  is 
computed  by  measuring  the  length  of  time  the  electric  field 
strength  is  above  a  threshold  set  by  the  user,  in  this  case 
.112  mV,  the  detection  threshold  used  in  the  NLDN  sensors. 


Signal  Strength  (LLP  Units) 


Fig.  5.5.  Signal  pulsewidth  vs.  range  for  the  30  kA  return 
stroke  modeled  in  Fig.  5.2.  Both  pre-  and  post-upgrade 
criteria,  10  ^is  and  7  |is,  are  depicted.  Note  that  the  range 
at  which  the  pulsewidth  drops  below  the  detection  criterion 
increases  when  we  decrease  the  pulsewidth  criterion. 
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Fig.  5.4  shows  that  the  peak  signal  strength  exhibits 
the  r~^  dependence  expected  from  the  radiation  field 
approximation  and  neglecting  attenuation.  The  RNSS  is 
constant  with  range,  also  as  expected  from  the  radiation 
field  approximation. 

Fig.  5.5,  the  pulsewidth  vs.  range  plot,  includes 
dotted  and  dashed  lines  at  10  p,s  and  7  jis,  respectively, 
corresponding  to  the  pulsewidth  detection  criteria  used  in 
the  NLND  DFs  before  and  after  the  upgrade  (refer  to  Chapter 
II) .  At  large  ranges  where  the  pulsewidth  is  less  than  the 
detection  criteria,  the  modeled  30  kA  return  stroke  would 
not  be  detected  by  the  model  DF.  The  range  where  the 
pulsewidth  equals  the  pulsewidth  detection  criteria  is  the 
theoretical  maximum  detection  range  of  the  DF  for  a  flash 
having  that  peak  current.  Fig.  5.5  shows  that  decreasing 
the  pulsewidth  detection  criteria  from  10  |j,s  to  7  |j,s 
increases  the  theoretical  maximum  detection  range  of  the  30 
kA  flash. 

f.  Repetitive  model  runs 

The  results  presented  thus  far  are  only  for  a  return 
stroke  with  peak  current  of  30  kA.  If  we  run  the  model 
repetitively  for  return  strokes  of  different  peak  currents 
we  can  derive  a  maximum  detection  range  vs.  peak  current 
relationship.  Fig.  5.6  depicts  this  relationship 
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Model  Detection  Range  vs.  Peak  Current 


Fig.  5.6.  Detection  range  vs.  peak  current  for  modeled 
return  strokes  having  waveforms  identical  to  Fig.  5.2  using 
detection  pulsewidth  criteria  of  10  ^s  and  7ns.  Note  that 
detection  range  for  a  flash  of  any  given  peak  current 
increases  when  we  decrease  the  pulsewidth  criterion. 
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for  return  strokes  of  the  same  waveform  and  peak  currents  of 

2,  10,  30,  50,  70,  90,  110,  and  130  kA  and  both  pulsewidth 
criteria. 

Changing  the  pulsewidth  criterion  has  no  effect  on  the 
RNSS  vs.  peak  current  relationship.  However,  we  see  in  Fig. 
5.6  that  decreasing  the  pulsewidth  criterion  increases  the 
theoretical  maximum  detection  range  of  a  DF  at  all  peak 
currents.  We  next  model  the  detection  of  return  strokes  by 
a  DF  to  investigate  the  effects  of  the  increase  in  detection 
range  on  the  detected  flash  peak  current  distribution. 

3.  Stochastic  return  stroke  detection  model 

Appendix  D  contains  the  MathCAD  6.0+  file  used  to  model 
return  stroke  detection.  Using  the  detection  range  vs.  peak 
current  relationships  derived  above  for  each  pulsewidth 
criterion,  the  detection  model  computes  peak  current 
distributions  of  detected  flashes.  This  is  accomplished  by 
generating  a  large  number  of  model  return  strokes  with 
random  peak  currents  and  random  locations  relative  to  a 
network  of  DFs  and  determining  how  many  are  detected  using 
the  pre-  (10  ns)  and  post-upgrade  (7  ^s)  pulsewidth 


criteria. 


Fig.  5.7.  Detection  model  randomly  generated  return  stroke 
locations  (only  1/10  of  total  shown  for  clarity) . 
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a.  Geometry  and  peak  current  distribution  of  the 
stochastic  model 

The  model  allows  the  user  to  set  the  number  and  x-y 
positions  of  multiple  model  DFs.  For  the  following 
discussion,  two  DFs  are  used,  located  along  the  x-axis  at 
x=+/-125  km.  This  250  km  spacing  between  the  two  model  DFs 
corresponds  to  the  approximate  average  spacing  of  DFs  in  the 
NLDN. 

10^  negative  return  strokes  are  generated  with 
locations  normally  distributed  about  the  origin  with  x  and  y 
standard  deviations  of  500  km.  Thus  we  expect  roughly  95% 
of  the  flashes  to  be  within  1000  km  of  the  origin.  Fig.  5.7 
depicts  one  tenth  of  the  randomly  generated  flash  locations 
to  provide  a  sense  of  the  flash  locations. 

The  10^  flashes  also  have  randomly  generated  peak 
currents.  The  peak  currents  were  generated  using  both  a 
distribution  and  a  uniform  distribution.  The  x^  peak 
current  distribution  is  probably  more  realistic  but  has  very 
small  sample  sizes  near  the  tails  of  the  distribution. 

Since  the  analysis  from  Chapter  IV  dictates  that  we  focus 
attention  on  weak  flashes,  the  small  sample  size  of  the  x^ 
distribution  near  0  kA  is  a  drawback.  While  not  as 
realistic,  the  uniform  distribution  provides  a  large  sample 
size  at  all  peak  currents  of  interest  and  eliminates  any 
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bias  the  initial  peak  current  distribution  has  on  the 
detected  flash  peak  current  distribution. 

b.  Detection  criteria 

We  consider  a  flash  to  be  detected  by  an  individual  DF 
if  it  is  within  the  maximum  detection  range  for  its  peak 
current/  as  determined  by  the  relationship  computed  earlier. 
In  the  NLDN,  however/  detection  by  multiple  DFs  is  required 
in  order  to  determine  the  location  of  a  flash.  We  therefore 
require  a  model  flash  to  be  within  detection  range  of  both 
DFs  to  be  detected. 

c.  Results 

Fig.  5.8  depicts  the  peak  current  distribution  of  all 
flashes  (solid  line) ,  and  the  peak  current  distribution  of 
flashes  detected  with  both  the  10  ps  (dashed)  and  7 
(dotted)  pulsewidth  criteria.  At  large  peak  currents/  above 
approximately  70  kA/  nearly  all  flashes  are  detected  with 
either  criterion.  At  lower  peak  currents,  however, 
significantly  more  flashes  are  detected  using  the  7  pis 
pulsewidth  criterion.  This  increase  is  maximized  around  20 
kA,  and  at  the  other  end  of  the  peak  current  spectrum,  near 
0  kA,  no  flashes  are  detected  using  either  pulsewidth 
criterion. 

The  top  half  of  Fig.  5.9  depicts  the  detection 
efficiency,  defined  as  the  fraction  of  flashes  detected,  for 


Model  Flash  Count  vs  Peak  Current 
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Fig.  5.8.  Modeled  return  stroke  detected  flash  count  vs 
peak  current  and  detected  flash  counts  at  10  jis  and  7  jis 
pulsewidth  criteria. 
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Fig.  5.9.  Model  detection  efficiency  vs.  peak  current  at  10 
|is  and  7  jxs  pulsewidth  criteria  (top)  and  model  post-upgrade 
flash  count  vs.  peak  current  normalized  to  pre-upgrade  flash 
count  (bottom) . 
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each  pulsewidth  criterion/  and  the  difference  in  detection 
efficiencies.  At  high  peak  currents,  both  pulsewidth 
criteria  result  in  near  100%  detection  efficiency  due  to  the 
very  large  maximum  detection  range  for  those  flashes 
compared  to  the  DF  spacing.  As  peak  current  decreases,  the 
detection  efficiency  falls  off  more  rapidly  for  the  10  ^is 
pulsewidth  criteria.  At  0  kA,  both  criteria  result  in  0% 
detection  efficiency.  The  pre-  to  post-upgrade  detection 
efficiency  increase  is  maximized  at  22  kA. 

The  bottom  half  of  Fig.  5.9  depicts  the  post-upgrade 
flash  count  vs.  peak  current  normalized  to  the  pre-upgrade 
value.  Comparing  the  bottom  half  of  Fig.  5.9  with  Figs. 

4.18  and  4.20  reveals  some  qualitative  similarities.  Both 
the  actual  and  the  model  post-upgrade  flash  count  ratios  are 
approximately  1  for  large  peak  currents.  This  ratio 
increases  to  a  peak  of  many  times  the  pre-upgrade  flash 
count  below  10  kA,  then  decreases  toward  0  as  peak  current 
decreases  toward  0  kA. 

In  the  US  NLDN  data,  the  post-upgrade  negative  flash 
count  peaks  at  7  times  the  pre-upgrade  count  at  6  kA  and  is 
approximately  equal  to  the  pre-upgrade  count  for  peak 
currents  above  30  kA.  The  modeled  post-upgrade  negative 
flash  count  peaks  at  18  times  the  pre-upgrade  count  at  8  kA 
and  is  approximately  equal  to  the  pre-upgrade  count  for  peak 
currents  above  50  kA.  The  mean  peak  current  of  the  model- 
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detected  return  strokes  decreases  from  65.3  kA  before  the 
hypothetical  upgrade  to  57.4  kA  after. 

Although  the  results  are  not  presented  here,  the  use  of 
the  peak  current  distribution  does  not  qualitatively 
affect  the  results. 


96 


CHAPTER  VI 

DISCUSSION  AND  CONCLUSIONS 
1.  Statistical  significance  of  data 

If  we  assume  that  the  yearly  mean  peak  current  and 
flash  counts  are  normally  distributed  about  pre-upgrade  mean 
values,  then  approximately  2.3%  of  the  yearly  mean  values 
will  lie  more  than  two  standard  deviations  above  and  2.3% 
will  lie  more  than  two  standard  deviations  below  the  pre¬ 
upgrade  mean  (Ott  1993) .  A  two  standard  deviation  increase 
or  decrease  in  mean  peak  current  or  flash  count  after  the 
upgrade  therefore  allows  us  to  conclude  with  97.7%  certainty 
that  the  post-upgrade  observed  value  is  not  caused  by 
natural  variance  in  the  data.  For  the  purposes  of  this 
discussion,  a  two  standard  deviation  increase  or  decrease 
will  define  a  significant,  artificial  change  in  the  data. 

As  discussed  in  Chapter  IV,  the  US  post-upgrade 
negative  flash  count  lies  within  two  standard  deviations  of 
the  pre-upgrade  negative  mean  peak  current,  so  we  cannot 
conclude  with  much  certainty  that  the  observed  increase  in 
negative  flash  count  has  other  than  natural  causes.  By 
contrast,  the  US  post-upgrade  positive  flash  count  is 
approximately  six  standard  deviations  higher  than  the  pre¬ 
upgrade  mean,  leading  us  to  conclude  that  the  observed 
positive  flash  count  is  due  to  factors  other  than  natural 


variability. 
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The  US  post-upgrade  negative  mean  peak  current  is 
approximately  three  standard  deviations  lower  than  the  pre¬ 
upgrade  mean,  leading  us  to  conclude  that  the  observed 
change  is  not  a  natural  phenomenon.  Finally,  the  US  post¬ 
upgrade  positive  mean  peak  current  is  approximately  five 
standard  deviations  below  the  pre-upgrade  value,  also 
leading  us  to  conclude  that  something  other  than  natural 
variability  is  responsible  for  the  observed  change. 

In  siommary,  then,  we  conclude  that  since  the  upgrade 
the  NLDN  is  detecting  significantly  more  positive  flashes 
and  that  the  mean  peak  current  of  detected  negative  and 
positive  flashes  is  significantly  lower. 

2.  Detection  modeling  results 

The  return  stroke  radiation  model  presented  in  Chapter 
V  shows  that  decreasing  the  pulsewidth  detection  criterion 
used  by  the  DFs  in  the  NLDN  effectively  increases  the 
maximum  detection  range  of  a  return  stroke  with  any  given 
peak  current. 

The  return  stroke  detection  model  further  shows  that  if 
the  maximum  detection  range  of  return  strokes  is  increased 
in  a  network  of  DFs  with  constant  spacing,  the  increase  in 
flash  count  occurs  primarily  for  weak  flashes.  Using  a 
random  sample  of  return  strokes  with  normally  distributed 
positions  and  uniformly  distributed  peak  currents,  post-  to 
pre-upgrade  flash  count  ratios  similar  to  those  found  in  the 
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archived  NLDN  data  result.  Additionally/  the  mean  peak 
current  of  detected  flashes  decreases. 

3.  Possible  explanation  of  peak  current  decrease  and  flash 
count  increase 

Based  on  the  combined  results  of  Chapters  IV  and  V,  a 
possible  explanation  of  the  post-upgrade  mean  peak  current 
decrease  is  the  increased  detection  efficiency  of  weak 
flashes  in  the  NLDN,  caused  by  the  decrease  in  DF  detection 
pulsewidth  criterion  from  10  |j,s  to  7  )is  during  the  upgrade. 

Decreasing  the  pulsewidth  detection  criterion 
effectively  increases  the  detection  range  of  a  given  flash. 
The  detection  range  of  strong  flashes  exceeds  the  average 
spacing  between  DFs  using  either  criterion  and  therefore 
guarantees  that  multiple  DFs  will  detect  strong  flashes, 
regardless  of  pulsewidth  criterion  used.  Weak  flashes,  on 
the  other  hand,  may  have  detection  ranges  significantly  less 
than  DF  spacing  using  one  or  both  pulsewidth  criteria.  But 
by  decreasing  the  pulsewidth  criterion  and  thereby 
increasing  the  detection  range  vs.  peak  current  relation, 
more  weak  flashes  are  detectable  by  the  NLDN  after  the 
upgrade . 

So  we  observe  an  increase  in  weak  flash  count  and 
little  or  no  change  on  strong  flash  count.  The  net  effect 
of  the  increased  presence  of  weak  flashes  in  the  NLDN  data 
is  to  decrease  the  mean  peak  current  of  detected  flashes. 
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Decreasing  the  pulsewidth  detection  criteria  is  not  the 
only  sensor  change  that  would  produce  this  result.  In  fact, 
any  change  which  increases  DF  detection  range  for  a  given 
flash  would  have  the  same  effect.  Since  all  of  the  changes 
incorporated  into  the  DFs  during  the  upgrade  were  intended 
to  improve  detection  efficiency/  they  all  had  the 
complementary  effect  of  increasing  the  percentage  of  weak 
flashes  detected. 

Note  that  this  explanation  does  not  account  for  the 
larger  increase  in  weak  positive  flash  count  than  in  weak 
negative  flash  count.  A  possible  clue  to  the  difference  can 
be  found  in  the  nature  of  cloud  flashes.  Since  cloud 
flashes  move  charge  between  cloud  charge  regions,  they 
typically  move  negative  charge  upward  or  positive  charge 
downward  (as  a  consequence  of  the  double-dipole  cloud  charge 
structure) ,  Both  processes  result  in  emitted  radiation  from 
cloud  flashes  of  positive  polarity  (i.e.  positive  charge 
moving  downward) .  Cloud  flashes  tend  to  have  small  peak 
currents  (<10  kA) ,  and  have  pulsewidths  on  the  order  of  10 
|is  (Krider  et  al.  1979)  . 

So  cloud  flashes  with  unusually  long  (>10  ps) 
pulsewidths  may  meet  the  CG  flash  detection  criteria  and  be 
mistaken  for  CG  flashes.  Decreasing  the  DF  pulsewidth 
detection  criteria  from  10  ^s  to  7  ^is  would  allow  even  more 


cloud  flashes  to  be  accepted  by  the  CG  flash  detection 
criteria. 
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The  increased  percentage  of  positive  flashes  observed 
after  the  NLDN  upgrade  may  therefore  be  a  consequence  of 
increased  contamination  of  CG  flash  data  by  cloud  flashes 
caused  by  decreasing  the  pulsewidth  detection  criterion. 
Increasing  sensor  gain  my  also  increase  detection  of  the 
weak  cloud  flash  radiation  signal  that  previously  went 
largely  undetected. 

If  we  look  closely  again  at  Fig.  4.19,  we  see  that  the 
positive  peak  current  distributions  for  pre-upgrade  years 
are  bimodal.  There  is  a  large  peak  between  10  and  15  kA  and 
a  much  more  subtle  peak  between  40  and  45  kA.  The  two 
separate  peaks  in  the  distribution  hint  at  two  separate 
processes  taking  place,  i.e.  detection  of  both  cloud  flashes 
and  positive  CG  flashes.  Fig.  4.19  shows  that  the  increase 
in  positive  flash  count  is  largest  near  the  lower  peak. . 

So  in  addition  to  decreasing  the  mean  peak  current  of 
detected  positive  and  negative  flashes,  the  NLDN  upgrade  may 
have  had  the  unwanted  effect  of  increasing  the  contamination 
of  positive  CG  flash  data  by  cloud  flashes. 

Another  possibility  is  that  the  NLDN  is  now  able  to 
detect  a  new  type  of  positive  CG  lightning  waveform  (Brook 


1996)  . 
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4 .  Recommendations 

These  results  suggest  that  caution  must  be  used  in 
comparing  lightning  data  from  before  and  after  the  1994 
network  upgrade.  As  Ciommins  et  al.  (1995)  point  out, 
"variations  in  [detection  efficiency]  can  seriously  affect 
the  measured  peak  current  distribution,  as  well  as  the 
measured  percent  of  events  that  have  a  positive  polarity" . 
With  this  in  mind,  it  may  be  beneficial  to  attempt  to 
develop  a  correction  factor  to  convert  pre-upgrade  flash 
distributions  to  equivalent  post-upgrade  values  in  order  to 
eliminate  the  effects  of  the  upgrade  in  comparing  annual 
data  sets. 

A  more  extensive  detection  model  that  incorporates  the 
actual  DF  positions  would  be  useful  in  determining  the  post¬ 
upgrade  effects  on  detected  flash  count  distributions  using 
the  actual  NLDN  DF  locations.  This  would  enable  an  estimate 
of  the  actual  NLDN  detection  efficiency  before  and  after  the 
upgrade,  and  provide  a  basis  for  such  a  conversion. 

Finally,  further  DF  calibration  studies  are  needed. 

The  existing  calibrations  by  Orville  (1991)  and  Idone  et  al. 
(1993)  consider  only  negative  triggered  return  strokes  with 
peak  currents  less  than  60  kA  detected  by  a  handful  of  DFs 
in  a  limited  geographical  area.  This  allows  the 
conductivity  of  Florida  soil  and  possible  azimuthal  biases 
of  the  DFs  used  in  the  calibration  to  skew  the  results.  A 
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new  post-upgrade  calibration  including  a  wide  geographic 
range  within  the  NLDN/  a  wider  peak  current  range,  and  both 
negative  and  positive  flashes,  would  improve  the  accuracy  of 
NLDN  peak  current  measurements. 
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APPENDIX  A 

ORVILLE  AND  IDONE  ET  AL.  PEAK  CURRENT  CALIBRATIONS 

This  appendix  contains  the  MathCAD  6.0+  files  used  to 
duplicate  the  Orville  1991  and  Idone  et  al.  1993  DF 
calibrations,  correct  the  Orville  1991  calibration  for  an 
error  in  the  position  of  one  of  the  DFs,  and  recalculate 
each  with  the  intercept  of  the  best-fit  line  constrained  to 


0  kA. 
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Verification  of  Calibration  Paper  by  Orville  (1991)  25  June  96 

Constants:  ORIGINsl  kA  lO^  amp  LLP:=l  tesla  (LLP  units  are  proportional  to  teslas) 


DF  ranges  from  trigger  site: 


range  := 


117.9  i 

LA 

145.1 

FP 

197.2 

•km  J' 

259.1 

Im 

379.4 

Sa 

426.8 

Qu 

(note:  this  range  is  incorrect) 


Read  in  data: 

French  measured  peak  current  (kA): 

OF  signal  amplitudes  (LLP  units): 
Normalized  DF  amplitudes  to  100  km: 


data  :=  READPRN(dataorv) 

nflashes  :=  rows(data)  i  :=  1  ..nflashes 


‘french.  =  data.  ,  kA 


ndfs  :=  cols(data)  -  1 


j  :=  1 ..  ndfs 


dfamp.  .  :=  data.  .  -LLP 

M.J  >0  + 1 


SI  range: 

dfamp.  . - 

■’J  lOO  kn 


normamp.  := 


^(dfampij.^^O) 
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Linear  regression: 


xbar  := 


^  'normamp. 
i 

nflashes 


ybar  := 


french. 

i 

nfrashes 


Sxx  :=  ^normampj  -  xbar^^ 
i 

Sxy  :=  ^  ^normampj  i  frgncjj.j  -  nflashes- xbar- ybar 
i 

syy 


Slope,  intercept,  and  best-fit  line: 


Correlation  coefficient  and 
standard  deviation: 


Sxx 


a  :=  ybar  -  xbar-b 


yhat(x)  :=  a  -r  b-x 


r  := 


Sxy 


Sxx-Syy 


S  := 


french; 


nflashes  -  2 


Measured  peak  current  vs.  DF  signal  strength  (normalized  to  100  km) 


Best-fit  intercept:  a  =  2.275  "kA 

kA 


slope:  b  =0.186' 


LLP 


correlation  coefficient:  r = 0.893 

standard  deviation:  s  =6.132*kA 


measured  peak  current  (kA) 


Consider  driving  the  intercept  to  0: 


yhat(x)  :=  b  x 


Measured  peak  current  vs.  DF  signal  strength  (normalized  to  100  km) 


DF  normalized  signal  strength  (LLP  units 


Ill 


Verification  of  Calibration  Paper  by  Idone  (1 993)  25  June  96 

Constants:  ORIGIN  kA  :=  lO^  amp  LLP  :=!  •  tesla  (LLP  units  are  proportional  to  teslas) 


DF  ranges  from  trigger  site: 


range 


117.9 

LA 

145.1 

FP 

220 

•km 

259.1 

Im 

379.4 

Sa 

426.8 

Qu 

Jk  (note:  this  range  is  corrected) 


Read  in  data: 

French  measured  peak  current  (kA): 

DF  signal  amplitudes  (LLP  units): 
Normalized  DF  amplitudes  to  1 00  km: 


data  READPRN(dataido) 

nflashes  :=  rows(data)  i  ;=  1 ..  nflashes 

*  french.  ”  data.  j  kA 

ndfs  :=  cols(data)  -  1  j  1 ..  ndfs 

dfamp.  .  :=  data.  .  ,  LLP 

i.j  +  1 


normamp. 


rangCj 
100  km, 


^(dfamp;  .^0) 


j 
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Slope,  intercept,  and  best-fit  line; 


Correlation  coefficient  and 
standard  deviation; 


Measured  peak  current  vs.  DF  signal  strength  (normalized  to  100  km) 


DF  normalized  signal  strength  (LLP  units 


Best-fit  intercept:  a=4.208-kA  correlation  coefficient:  r  =0.881 


slope:  b  =0.171  •— — 

LLP 


standard  deviation:  s  =4.649  -kA 


Consider  driving  the  intercept  to  0: 


2j“°nnamPii  french, 
i 

irA 

K  —  n  ocw  • 

^  ^normampj^^ 

i 

LLP 

Sxy 

2  f  frenchj  -  yhat(nonnampj)j2 

Jsxx  Syy  s  := 

i 

nflashes  -  2 

yhat(x)  :=  b  x 


Measured  peak  current  vs.  DF  signal  sfrength  (normalized  to  1 00  km) 


s=  4.984 -kA 


Correction  to  Calibration  Paper  by  Orville  (1991)  25  June  96 

Constants:  ORIGIN=l  kA  ;=  lO^amp  LLP- 1-tesla  (LLP  units  are  proportional  to  teslas) 

LA 
FP 

Jk  (note:  this  range  is  corrected) 

Im 
Sa 

Qu 

Read  in  data:  data  ;=  READPRN((lataorv) 

nflashes  :=  rows(data)  i  :=  1 ..  nllashes 
French  measured  peak  current  (kA):  i 

ndfs  -  cols(data)  -  1  j  :=  1 ..  ndfs 

DF  signal  amplitudes  (LLP  units):  dfamp.  .  :=  data.  .  ^  j-LLP 

Normalized  OF  amplitudes  to  100  km: 

S 

normamp.  := 

^(dfampj.^O) 

j 
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Linear  regression; 


xbar  := 


nonnamp. 


nflashes 


ybar  := 


‘  french. 


nflashes 


Sxx  :=  ^  ^normampj  -  xbar^^ 
i 

Sxy  ;=  ^  ^normampj  i  g-endj  j  -  nflashes  xbarybar 
i 


Slope,  intercept,  and  best-fit  line: 


a  :=  ybar  -  xbar  b 


y(x)  :=  a  +  b  x 


Correlation  coefficient  and 
standard  deviation: 


X](‘ french; 


JSxxSyy 


nflashes  -  2 


Measured  peak  current  vs.  DF  signal  strength  (normalized  to  100  km) 


)0  150  200  2: 

DF  normalized  signal  strength  (LLP  units 


300  350 


Best-fit  intercept:  a  =  2.684  'kA 


correlation  coefficient:  r  =0.887 


slope:  b=0.178* 


standard  deviation:  s  =  6.304  -kA 


measured  peak  current  (kA) 


Consider  driving  the  intercept  to  0: 


b:= 


X!“onnamPiifi«ichj 

i 


^normampj^ 


b  =0.195 


kA 


LLP 


yhat(x)  :=  b  x 


r  := 


Sxy 


-JSxxSyy 


S  := 


E  french;  -  yhat^nonnamp;^^^ 


nflashes  -  2 


Measured  peak  current  vs.  DF  signal  strength  (normalized  to  1 00  km) 


s=  6.427 -kA 
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APPENDIX  B 

NLDN  DATA  ANALYSIS  SAMPLE 

This  appendix  contains  an  example  of  the  MathCAD  6.0+ 
files  used  to  perform  the  statistical  analysis  presented  in 
Chapter  IV.  A  separate  file  was  used  for  positive  and 
negative  flashes  for  Kansas,  the  Ohio  Valley,  and  Florida, 
and  the  entire  continental  United  States.  FLASH  and  the 
FORTRAN  code  mentioned  in  Chapter  IV  were  used  to  obtain  the 
flash  counts,  mean  peak  currents,  and  peak  current 
distributions  processed  by  these  MathCAD  files. 
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Statistical  Summary  of  U.S.  Negative  Flashes,  1 989-95  2e  Jui  96 

Read  in  data  and  compute  totals,  means,  standard  deviations,  and  distributions 

ORIGIN^  1 

n  :=  READPRN(usnneg)  number  of  flashes  by  month  and  year 

years  :=  cols(n)  y  1 ..  years 

months  :=  rows(n)  m  :=  1 ..  months 
yearsbefore  :=  years  -  2 

ntot  :=  >  n  total  annual  flash  count 

y  «n,y 

m 

yearsbefore 

^  y=l  pre-upgrade  annual  mean  flash  count 

annual  =  yearsbefore 


nsd 


annual 


yearsbefore 


E  annual)^ 


y=l 


yearsbefore 


pre-upgrade  annual  flash 
count  standard  deviation 


month 

m 


yearsbefore 

E  n 

Z_i  “-y 

y=i 

yearsbefore 


pre-upgrade  monthly  mean  flash  count 


month 

m 


yearsbefore 

y=  1 _ 

yearsbefore 


pre-upgrade  monthly  flash 
count  standard  deviation 
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i  :=  READPRN(usin^) 


mean  peak  current  by  month  and  year 


m,y 


annual 

y 


annual  mean  peak  current  (weighted  by  flash  count) 


annual  peak  current  standard 
deviation  (weighted  by  flash  count) 


month  = 

m 


yearsbefore 


y=l 


yearsbefore 

Z_j  ®.y 

y=l 


pre-upgrade  monthly  mean  peak 
current  (weighted  by  flash  count) 


month  ■ 
m 


yearsbefore 

y~!  °m,  y'  y  ~  month^ j 


y=l 


yearsbefore 


y  n  -  1 
Z_i  “.y 

y=l 


pre-upgrade  monthly  peak 
current  standard  deviation 
(weighted  by  flash  count) 
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dist  :=  READPRN(usdistn) 
bins  :=  4  j  :=  1 ..  bins 

“winter^  = 

°spring._y'=‘^®^j  +  4.y 
“summer  *®*j  +  *,y 

“  autunuij  y  ”  j  + 12,  y 

n*.*  :=  dist.  ,, 

j.y  J  + 


P  winter. 

J 


“winter. 

J.y 


^“winter. 

”  J.J 


‘  spring 


P  snri 


J.y 


sprmg. 

Xi”sP““gj,j 

j 


P  summer. 

J.y 


summer. 

J.y 


E” 


summer. 

J.y 


“autumn. 


P  autumn. 


J.y 


J.y  >  n 


/  .“autumn. 

—  J.) 

J 


peak  current  distribution  by  season 


matrix  index 


seasonal  flash  count  distributions 

i  Mq 

1  <20  kA 

2  20-30  kA 

3  30-40  kA 

4  >40  kA 


total  annual  flash  count  distribution 


seasonal  flash  count  distributions 
expressed  as  a  fraction  of  the  annual  total 


number  of  flashe«  number  of  flashes 


Graphical  Results 


U.S.  Annual  Total  Negative  Flash  Count 


Note:  dashed 
line  is  2 
standard 
deviations 
above  and 
below  the 
pre-upgrade 
annual  mean 


1988 


1989 


1990 


1991 


1992 


1993 


1994 


1995 


1996 


year 


U.S.  Monthly  Negative  Flash  Count 


S'lO”  X 


6-10®  -f- 


4*10  + 


2*10''  + 


T  / 
/ 


.-y 

y/ 


:--r- 


--•'f 


dotted  =  1994 
dashed  =  1995 


Note:  error 
bars  are  two 
standard 
deviations 
above  and 
below  the 
pre-upgrade 
monthly  mean 


•N 

•v 


10 


11 


12 


8 


9 


mean  peak  current  (kA)  mean  peak  current  (kA) 
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l 


U.S.  Annual  Mean  Negative  Peak  Current 


Note: 
error  bars 
are  two 
standard 
deviations 
above 
and  below 
the  mean 


U.S.  Monthly  Negative  Mean  Peak  Current 


80  T 


60  + 


40 


20  + 


i-i-i 


...I 


dotted  =  1994 
dashed  =  1995 


•I- . 


Note:  error 
bars  are  two 
standard 
deviations 
above  and 
below  the 
pre-upgrade 
monthly  mean 


10 


11 


12 


month 


fraction  of  flmhes  fraction  of  flashes 
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U.S.  Spring  Negative  Peak 
Current  Distribution 


U.S.  Summer  Negative  Peak 
Current  Distribution 


U.S.  Autumn  Negative 
Peak  Current  Distribution 


U.S.  Winter  Negative  Peak 
Current  Distribution 
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Supplemental  stats  (11/25/96): 


yearsbefore 

Y,  *‘^annuay“‘oty 

y=  I 

ibarpre  := -  isdpre  := 

yearsbefore 

S  ““'y 

y-l 

^  ntoty-nbar  annual  annual^  “ ‘^'Pre 

^  annual 


yearsbefore 

Y  [“toy  (ibar  ^nnual^  “  ibaiprej^ 

y  =  i _ 

yearsbefore 

ntoty  -  1 

y=l 


1.3*  lo’ 


ntot  = 


1.532*  lo’ 
1.623*  lo’ 
1.558*10’ 
2.313*10’ 
2.302*  lo’ 
2.061  *10’ 


ntot  =1.269*10* 
y 


y 

°barannual  =  b665*10’ 

annual  =3415772.271 
annual  =6831544.542 


-1.069 


Zn  = 


-0.389 

-0.124 

-0.314 

1.895 

1.865 

1.159 


ibar 


annual 


37.957 

2.425 

41.574 

2.853 

37.657 

2.446 

35.821 

annual 

2.316 

35.626 

3.59 

34.259 

1.727 

30.249 

2.267 

0.205 

1.893 

ibarpre  =37.517 

0.065 

Zi  = 

-0.791 

isdpre  =2.144 

-0.882 

2isdpre  =4.287 

-1.52 

-3.391 
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Finite  Difference  Lightning  Return  Stroke  Field  Model  02  Oct  96 

Units  and  Constants 

kA  :=  10^-amp  ps  :=  10‘^  sec 

Electromagnetic 
constants 

LLP  Units 

Return  stroke  peak  current  and  wavefront 
velocity 

Exponentially  decaying  current  after  initial 
‘mod(^>‘)  =  [‘peak  e  J  (h<v  t)  breakdown  peak 


Return  Stroke  Current  Model 

ipeak  :=  30  kA  v  :=  10*  — 

sec 


8  m 

c:=:310* - 


.  .  7  newton 

H0:=  4-7c10^ - - 


60=  8.85- 10' — 
newton- 


LLP  3  =  1-510 


10  weber 


LLP  g  c-LLP  3 


nh  ;=  40  Channel  top  and  bottom  heights  and 

number  of  steps  for  integration 

bj  ;=  h  mjjj  +  j-Ah  Channel  step  size:  Ah  =  50*m 

Initial  return  stroke  channel  charge  density 


Channel  Current  and  Charge  Density 

nt:=1000 


Time  limit  and  number  of  time  steps  for 
channel  current  density  computation 


^;=k-At  Time  step  size:  At =0.05 -ps 

Courant  condition:  v- At  =  5-m 

Channel  current  (constrained  to  0  at  top  of 
channel) 
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^j.o 


=  Poft) 


Initial  channel  charge  density 


k:=  l..nt 


n.k  - 


\k-i  Vk-i 


^j.k-  1 

Pj.k-  1 

Pj.k-  1  ■ 


Ah 


/i 


‘nb,k-  1  *nh-  l,k-  1 


Ah 


fj  +  l,k-l  j-l,k-l 


At  if  j«0 

At  if  j=nh 


2Ah 


At  otherwise 


Channel  charge  density  (determined  by 
continuity  from  the  vertical  derivative  of 
channel  current) 


k  :=  300  tj,  =  15*^s 


Channel  Current  Channel  Charge  Density 


Vector  and  Scalar  Potentials 


R(r 


,h)  :=  Jr^  + 1 


kret(r,t,h)  := 


0  if  t- 


floor 


R(r,h) 


\  R(r.h)\ 


At 


<0ns 


otherwise 


Slant  range  expression 

Retarded  time  index  (accounts  for  time 
required  for  radiation  to  travel  from  a  spot 
on  the  channel  to  the  observer) 


Vertical  component  of  the 
vector  potential  (factor  of  2 
accounts  ofr  image  current) 
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•Kr.z.t) 


nh 

E 

Fj,kret^r,t,h.  -  zj 

/Ah' 

Rfr.h.  -  z) 

Rfr,h.  ,  -  z) 

j  =  i 

Fj,kret^r,t,b.  +  z 

V  J- 1  J 

\ 

/ 

+ 

"Fj-  l,kret^r,t,b._  j  +  zj 

.1 

^Ah 

R(r.k.,z) 

R(r.h._j  +  z) 

~2 

Azimuthal  Magnetic  and  Vertical  Electric  Fields 


Scalar 

potential 

(sum  within 

summation 

accounts  for 

image 

charge) 


At  :=  .2  km 
Az  :=  1-km 
At  :=  .5jis 


Radial  distance  step  for  spatial  derivative  of  vector  potential 
Vertical  distance  step  for  spatial  derivative  of  scalar  potential 
Time  step  for  temporal  derivative  of  vector  potential 
Courant  condition:  c-  At = 0. 1 5  -km 


B<|)(r.t)  :=  - 


/Az(r  + Ar,t)  -  Az(r-  Ar,t)\ 

\  ^A^  I 


Azimuthal  component 
of  the  magnetic  field 


Ez(r,t)  :=  - 


|(i>(r.Az,t)  -  4)(r,-Az,t)j 


/Az(r,t  At)  -  Az(r,t  -  At)\ 

\  2^At  I 


Detection  Criteria 


Vertical  component  of 
the  electric  field 


time(r,  thresh) 


time<-0•^ls 

r 

t« — 
c 


while  t<(-  +  40- ps 
c 


(time^-time  +  Ips) 
t«-t+  Ips 


time 


if  (Ez(r,t)>thresh) 


Amount  of  time  detected 
signal  is  above  a  given  signal 
strength  threshold  at  a 
particular  range 


peak(r)  := 


peak^-O 


weber 


m 


r 

t< — 
c 


while  t<|-  -t-  40- ps 
c 


(peak«-B(ti(r,t))  if  (B4i(r,t)>peak) 
t«-t  +  Ips 


peak 


Peak  (magnetic)  signal  strength 
of  detected  radiation 


r 

100  km 


mss(r)  :=  peak(r)- 


Range-normalized  signal  strength  (peak 
magnetic  field  strength  normalized  to 
100km) 
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APPENDIX  C 

FINITE  DIFFERENCE  FIELD  MODEL 

This  appendix  contains  the  MathCAD  6.0+  file  used  in 
the  return  stroke  radiation  model  discussed  in  Chapter  V. 

The  user  supplies  the  return  stroke  model  characteristics  as 
a  function  of  time  and  height  and  the  model  computes  the 
electromagnetic  fields  radiated  by  the  return  stroke,  the 
raw  and  normalized  signal  strength  vs.  range,  and  the  signal 
pulsewidth  vs.  range.  Multiple  runs  of  this  model  using 
return  strokes  of  various  peak  currents  produces  the 
detection  range  vs.  peak  current  relation  used  in  the 


detection  model. 
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APPENDIX  D 

STOCHASTIC  DETECTION  MODEL 

This  appendix  contains  the  MathCAD  6.0+  file  used  to 
model  return  stroke  detection  with  pulsewidth  detection 

criteria  of  10  |is  and  7  ^is.  It  generates  a  large  number  of 
randomly  located  return  strokes  with  random  peak  currents 
and  determines  the  flash  count  distributions  of  detected 
flashes  using  each  detection  criterion. 
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Return  Stroke  Radiation  Detection  Model  04  Oct  96 

Units  and  Constants 

kA  :=  lO^  amp  LLP  g  1.5- 10'  lo. ORIGIN^  1 


Data  from  Previous  Finite  Difference  Model  Runs 


peak  cun'ents 


ranges  and  range-normalized  signal  strengths 


ipeak  := 


0 

40.811 

64.096 

2 

252.189 

399.107 

10 

631.957 

1.001-10^ 

30 

50 

l.OOMO^ 

2.513-10^ 

70 

•kA  mglO  ;= 

1.586- 10^ 

•km  mg7  = 

2.513-10^ 

90 

1.586- 10^ 

3.982-10^ 

110 

2.513- 10^ 

3.982- 10^ 

130 

■J 

2.513-10^  J 

6.311-10^ 

mss  10 


8.496 

8.475 

42.24 

42.223 

126.638 

126.617 

211.029 

-LLPg  mss7  := 

210.994 

295.41 1 

295.392 

379.814 

379.774 

464.187 

464.169 

548.585 

548.549 

i  :=  1 ..  rows(ipeak) 


counting  index  for  data 


(Ion)  RNSS  (LLP  units) 


Graphical  Results 


Signal  Strength  vs  Peak  Current 


peak  current  (kA) 


Best-fit  parameters  for  lOus  threshold 

LLPfi 

slope(ipeak,msslO)  =4.219* - 

kA 


mtercept(ipeak,nisslO)  =0.053  "LLPg 


Best-fit  parameters  for  7us  threshold 

LLPg 

slope(ipeak,mss7)  =4.219* - 

kA 


intercept (ipeak,mss7)  =0.033  *LLP  g 


Best-fit  parameters  for  10us  threshold 

km 

slope(ipeak,mglO)  =20.062* — 

kA 

intercept(ipeak,mglO)  =31.667  *km 

Best-fit  parameters  for  7us  threshold 

Igu 

slope(ipeak,mg7)  =44.305  * — 

kA 

intercept(ipeak,mg7)  =-129.128  *km 


peak  current  (kA) 
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Model  Detection  Efficiency  Calculations 
Predefined  functions 

(ietectmglO(i)  :=  slope(ipeak,mglO)  i 

detectmg7(i)  :=  slope(ipeak,mg7)  i 
mg(xl  ,yl  ,x2,y2)  :=  J(xl  -  xlf  +  (yl  -  y2)^ 


Model  Parameters 


d&  := 


ndfs  :=  rows(d6c) 

nflashes  :=  10*  i  :=  1  ..nflashes 


X  :=  momi(nflashes, 0,500)  km 
y  :=  raorm(nflashes,0,500)km 


current.  :=  md(100)  kA 

Detection  Determination 
j  :=  1 ..  ndfs 


dist  j  :=mg(x;.yj,dfej,dfyj) 
detectlO.  := 


detect7.  := 


^dist  j  s  detectmg  10  ^current 

j 

^dist.  .  <detectmg7  ^current^^ 


current. 


current. 


Maximum  detection  range  as  a 
function  of  return  stroke  peak 
current  at  10  and  7us  thresholds 


Distance  formula 


DF  locations 

Number  of  dfs 
Number  of  flashes 

Flash  location  distribution 

Flash  current  distribution 


Distance  from  flash  to  each  df 


Determine  if  range  from  flash  to 
each  df  is  less  than  detection 
range  at  each  threshold— if  so, 
then  the  flash  is  detected 
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Flash  Detection  Histogram 
n  :=  1 ..  100 
ipeak^  :=  n  kA 

htot  :=  hist(ipeak,cunent) 
hlO  :=  hist(ipeak,dGtectlO) 
h7  :=  hist(ipeak, detect?) 


n:=  1..99 

hlO  h? 

efflO  := - eff7  = -— 1- 

htot  “  htot 


fiac  := 

n 


h? 


hlO 


if  hlO  5=0 


0  otherwise 


Peak  current  bins  for  histogram 

Total  number  of  flashes  vs  peak  current 

Number  of  flashes  detected  vs  peak  current  with 
1 0  and  7us  thresholds 


Detection  efficiency  expressed  as  percent  of 
flashes  detected  at  1 0us  and  7us  thresholds 


Number  of  flashes  detected  at  7us  threshold 
expressed  as  a  fraction  of  number  of  flashes 
detected  at  1 0us  threshold 


^~|hl0^ipeak^ 

meanlO  = 

n 

2h7„ipeak„ 

mean?  := 


Mean  peak  current  of  detected 
flashes  at  1 0us  threshold 

meanlO  =65.3  -kA 


Mean  peak  current  of  detected 
flashes  at  7us  threshold 


n 


mean?  =57.381  ’kA 
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Graphical  Results 


Number  of  Detected  Flashes 


solid  =  total  number  of  flashes 
dash  =  detected  at  1 0us  threshold 
dot  -detected  at  7us  threshold 


Detection  Efficiency 


x=  10us  threshold 

+  =  7us  threshold 

dash  -  difference  of  1 0us  and 

7us  threshold  efficiencies 


Fractional  Increase  in  Detected  Flashes 


Number  of  flashes 
detected  atthe  7us 
threshold  expressed  as  a 
fraction  of  the  number  of 
flashes  detected  at  the 
1 0us  threshold 


peak  current  (kA) 
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10,20..  10^ 


DF  and  Flash  Locations 


DF  locations  (x's)  and 
locations  of  every  tenth 
randomly  generated  flash 


x(km) 
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